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Abstract 

Hadrons carrying open heavy flavor, i.e. single charm or bottom quarks, are among 
the key diagnostic tools available today for the hot and dense state of strongly interacting 
matter which is produced in collisions of heavy atomic nuclei at ultra-relativistic energies. 
First systematic heavy-flavor measurements in nucleus-nucleus collisions and the reference 
proton-proton system at Brookhaven National Laboratory’s (BNL) Relativistic Heavy Ion 
Collider (RHIC) have led to tantalizing results. These studies are now continued and 
extended at RHIC and at CERN’s Large Hadron Collider (LHC), where considerably 
higher collision energies are available. This review focuses on experimental results on 
open heavy-flavor observables at RHIC and the LHC published until July 2012. Yields of 
heavy-flavor hadrons and their decay products, their transverse momentum and rapidity 
distributions, as well as their azimuthal distributions with respect to the reaction plane in 
heavy-ion collisions are investigated. Various theoretical approaches are confronted with 
the data and implications for the properties of the hot and dense medium produced in 
ultra-relativistic heavy-ion collisions are discussed. 
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Introduction 


1.1 Strongly interacting matter in ultra-relativistic heavy-ion colli¬ 


sions 


Quantum chromodynamics (QCD) is the underlying theory of the strong interaction [I]. At 
sufficiently high temperature, T, or baryo-chemical potential, /in, QCD inspired model calcu¬ 
lations predict a phase transition of strongly interacting matter from a system of hadrons to 
a deconfined medium, dubbed a quark-gluon plasma (QGP), in which the relevant degrees of 



lattice QCD calculations become significantly more difficult. Some calculations predict that 
the phase transition remains of the cross over type for all values of /i B . Several others indicate 
that towards higher /i B the phase transition will eventually be of first order. Consequently, 
this would imply the existence of a critical point in the phase diagram of strongly interacting 
matter. At very high /j, B the presence of more exotic phases, e.g. color superconducting forms 

A schematic QCD phase diagram is shown 


14 


16 


of strongly interacting matter, is expected 
in Fig. [T[ 

Strongly interacting matter at extreme temperature and/or density is not only of fundamen¬ 
tal interest in the context of mapping out the QCD phase diagram, but it is also of astrophysical 
relevance. The early universe is believed to have spent its first few microseconds after the big 
bang in the state of a quark-gluon plasma at high temperature and zero baryo-chemical poten¬ 
tial. Hadronization took place when the universe expanded and cooled down below Tq. QGP 
matter at high density and moderate temperature might be present today inside of neutron 
stars 17 , or it might be created briefly during the core collapse of a supernova explosion [l8,19]. 


Collisions of heavy atomic nuclei at ultra-relativistic energies provide the unique opportunity 
to investigate experimentally the properties and the dynamics of hot and dense QCD matter 
in the laboratory, because only in such collisions the necessary temperature or density can be 
reached. 

First “circumstantial evidence” for the production of a quark-gluon plasma was reported 20 


from the fixed target heavy-ion program at the CERN Super Proton Synchrotron (SPS), where 
collisions of lead nuclei (Pb+Pb) were investigated at energies per nucleon-nucleon pair up 
to ^/snn = 17.3 GeV in the center of mass frame. Experiments with colliding gold nuclei 
(Au+Au) at ^/snn = 200 GeV at the dedicated BNL Relativistic Heavy Ion Collider (RHIC) 
have substantiated these findings and have led to a first quantitative characterization of the 
properties of the QGP 121— 24 , revealing some surprises. Here, only some of the most striking 


observations are briefly summarized. 

The initially produced fireball has such a high temperature and density that the partons, 
i.e. quarks and gluons, equilibrate on a time scale of less than 1 fm/c. Large pressure gradients 
in the system lead to a hydrodynamic evolution of the fireball which subsequently expands, 
cools down, and undergoes hadronization. When the system is cold and dilute enough, the 
hadrons freeze out chemically and thermally and stream to the detectors. It is important to 
note that the geometrical shape of the initial equilibrated partonic fireball is asymmetric for 
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Figure 1: Schematic phase diagram of strongly interacting matter for two massless quarks as a 
function of the temperature T and the baryo-chemical potential /jb . 


collisions with non-zero impact parameter, exhibiting an almond-like shape averaged over many 
collisions. Consequently, the pressure gradients driving the expansion are asymmetric as well. 
Thus the initial spatial anisotropy is translated into an azimuthal anisotropy in momentum 
space of the produced hadrons. Originally, this anisotropy was quantified based on a Fourier 
expansion of the momentum distribution 


d 3 N d 3 N ~ „ 

d 3 p p t d<j)dp t dy 


*r)] 


( 1 ) 


where </> is the azimuthal emission angle of a particle with respect to the orientation of the 
reaction plane in a given collision 25 . This plane is defined by the beam direction and 


the impact parameter vector connecting the centers of the colliding nuclei in the transverse 
plane. While the short axis of the initial almond-like shape of the fireball is aligned with the 
reaction plane, the long axis is perpendicular to the latter. The harmonic coefficients, v n , 
quantify the strength of the anisotropy. The second harmonic, v 2 , is also called the elliptic 
flow, and it is believed to be sensitive to the early dynamics in the collision. Recently, it was 
recognized that fluctuations of the initial state geometry are important as they can generate 
higher nth-order flow components 26-32 , in particular odd-n components, when measured 


with respect to the proper initial-state planes of symmetry <h n 33 instead of <Fr. For a given 
initial state configuration and equation of state of the expanding medium the laws of ideal 
or viscous relativistic hydrodynamics predict the momentum space azimuthal asymmetry. The 
observed good agreement between the predictions |40-42 and the measured v 2 of various hadron 
species imply that the produced medium expands collectively almost like a perfect fluid, i.e. 
a fluid with a small shear viscosity to entropy density ratio, 77 /s, close to the conjectured 43 


quantum lower bound of 77 /s = 1/47T. Precision measurements of iq or even higher orders of the 
azimuthal anisotropy (34]- 39] should improve the sensitivity to the viscosity of the medium |30|. 

This partonic fluid produced in Au+Au collisions at RHIC is so dense that quarks and 
gluons from initial partonic hard scattering processes (“jets”) suffer a large energy loss while 
they propagate through the hot and dense medium. This “jet quenching” gives rise to a strong 
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suppression of hadrons originating mainly from jet fragmentation at high transverse momentum, 
p t , relative to the expectation from a superposition of independent binary nucleon-nucleon 
collisions. To quantify this suppression the nuclear modification factor Raa is introduced as 


Raa — 


dN A A/dp t 


dN A A/dp t 


(Neon) x dNpp/dpi {Taa) x dcTpp/ dpi 


( 2 ) 


where dNAA/dpt is the differential invariant yield in nucleus-nucleus collisions and dN pp /dp t 
( d<Tpp/dp t ) is the corresponding differential invariant yield (cross section) in pp collisions. For 
a given centrality class the average number of binary collisions is denoted by (N co u) and (Taa) 
is the corresponding nuclear thickness function, which relates (N co n) with the inelastic nucleon- 
nucleon scattering cross section. (iV co n) and (Taa) are obtained via Glauber model calcula¬ 
tions 44,45 of the collision geometry taking into account the response of the detectors used 


for the centrality measurement. As expected for particles not participating in the strong in¬ 
teraction, the nuclear modification factor of photons from initial hard scattering processes, e.g. 
quark-gluon Compton scattering, is consistent with one (461. Furthermore, this confirms that 
the Glauber model based calculation of the collision geometry is under control and that, at 
least at mid-rapidity at RHIC, initial state modifications of the parton distribution functions 
in nuclear matter, e.g. shadowing, anti-shadowing, or saturation effects, do not play a sig¬ 
nificant role. For high p t hadrons from jet fragmentation the value Raa ~ 0.2 measured at 


mid-rapidity in Au+Au collisions at RHIC 121— 24,47 demonstrates the opaqueness of the pro¬ 


duced dense fireball. When confronted with these data, theoretical models indicate that initial 


energy densities beyond 10 GeV/fm 3 are reached in the produced fireball 48,49 , clearly in the 
QGP regime. 

The initial temperature of the fireball can not be measured directly. However, an enhanced 
yield of quasi-real virtual photon^] above known hadronic sources, was measured [50 in central 
Au+Au collisions at RHIC. Interpreting this excess as virtual thermal photon emission from 
the hot fireball, hydrodynamical models allow to infer the initial temperature T illit and ther- 
malization time To of the fireball. Values of T init ~ 300 — 600 MeV and To ~ 0.6 — 0.15 fm/c 
are consistent with the data [5l( . 

The yield of heavy quarkonia in nucleus-nucleus collisions, i.e. bound states of charm or 
bottom quarks and their antiquarks, was long thought of as an observable directly related to 
deconfinement in the fireball 1521. It was predicted that color charge screening in a deconfined 
QGP would prohibit charm quarks and antiquarks to form a bound state, the J/xp meson, 
similar to the process of Debye screening in an electromagnetic plasma (521. The resulting 
J/xp suppression was indeed observed for the first time in central nucleus-nucleus collisions 
at the CERN SPS 1531. However, other mechanisms unrelated to deconfinement, e.g. the 
absorption of quarkonia in cold nuclear matter or their break-up by interactions with other 
hadrons (co-movers) produced in the collision, could also describe the SPS data reasonably 
well |54--j56| rendering the interpretation inconclusive. Subsequently, it was realized that in a 
deconfined medium a new mechanism for quarkonium production might become important if a 
sufficiently large number of heavy quark antiquark pairs is produced initially. Quarkonia could 
be formed in this case via statistical hadronization at the phase boundary (57, 58| or, earlier, 
via coalescence in the QGP 59 . While in nucleus-nucleus collisions at the SPS the number 


of charm quark anti-quark pairs is much smaller than one, leaving no room for a significant 
contribution to the J/xp yield via regeneration processes, the situation is different at RHIC, 


Mow-mass e + e pairs with high p t (m e + e - < 0.3 GeV/c 2 , 1 < p t < 5 GeV/c) 
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where more than 10 charm quark anti-quark pairs are produced in a central Au+Au collision 
and, in particular, at the LHC, where more than 100 such pairs are generated in a central Pb+Pb 
collision. Under such conditions J /ijj suppression might actually turn into an enhancement via 
regeneration processes |60|. First evidence in this direction was observed at RHIC in the 
centrality and, in particular, the rapidity dependence of the nuclear modification factor of 3/ip 
mesons, R^/a 61 . R&\ reaches its maximum value at mid-rapidity in central Au+Au collisions. 
This is at variance with the expectation from suppression via color screening or from destruction 
through the interaction with co-moving hadrons, which both should reach their maximum near 
mid-rapidity where the energy density reaches is the largest. In the regeneration picture, 
however, this maximum of Raa a t mid-rapidity finds its natural explanation in the fact that 
in this region of phase space the density of charm quarks and anti-quarks is maximal. Still 
the situation is not conclusive yet, as other possible scenarios might explain the RHIC data. 
Most likely, the J/^ saga will only be settled by decisive measurements from the LHC, where 
rather dramatic effects due to charmonium regeneration are expected when the total charm 
production yield is large enough |62|. 

At the CERN LHC colliding beams of lead nuclei are investigated at unprecedented high 
energies (^/snn = 2.76 TeV up to now) since November 2010. The measurement of two-pion 
Bose-Einstein correlations via Hanbury-Brown Twiss interferometry has shown that the volume 
of the produced fireball is about twice as large as it was observed in Au+Au collisions at RHIC 
and has a longer life time as well [631. The pseudorapidity density of charged particles produced 
at mid-rapidity also increases by a factor close to two going from RHIC to the LHC |64], 
indicating a substantial increase of the initial energy density reached in Pb+Pb collisions at 
the LHC. Remarkably, for various hadron species the elliptic flow strength v 2 as a function 
of p t is the same within uncertainties at RHIC and at the LHC |65]-|67|, which demonstrates 
that the produced fireball still resembles the properties of an almost perfect fluid, albeit with 
a large initial energy density. At the LHC, high energy jets in the TeV range are copiously 
produced. Such hard probes for the medium are not available at RHIC and, consequently, the 
measurement of the nuclear modification factor can be extended to much higher p t in Pb+Pb 
collisions at the LHC compared to Au+Au collisions at RHIC. In central collisions, Raa reaches 
its minimum value, corresponding to maximum suppression of hadrons from jet fragmentation, 
for a hadron p t of about 6-7 GeV/c |68, 69 . Raa increases towards higher values of p t but 
stays significantly below one even for hadrons with a p t of 100 GeV/c |69], highlighting the 
enormous density of the produced medium. Also in the quarkonium sector intriguing new 
observations have been made already. The nuclear modification factor of 3/ip measured at 


forward rapidity |70| is significantly larger in central Pb+Pb collisions at the LHC compared to 
Au+Au collisions at RHIC, which might indicate regeneration of charmonia from a deconfined 
medium. In addition, the suppression of various T states, i.e. bottomonia, can be studied in 
a systematic way at the LHC. After first indications for T suppression have been reported at 
RHIC 71 , where the different spin states have not been resolved yet, first data from Pb+Pb 


collisions at the LHC show a suppression of the Y(1S) ground state in central collisions at low 
Pt (I?) 

For the coming years the most important scientific challenge for the field is the quantitative 
characterization of the state of QCD matter produced in ultra-relativistic heavy-ion collisions at 
RHIC and at the LHC. In this context the measurement of observables related to the production 
of hadrons carrying open heavy flavor can play a unique role as is argued in Section [H 
where the basic concepts and related theoretical approaches are briefly introduced. Results, 
relevant for the current discussion, originating from open heavy-flavor measurements in hadronic 
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Table 1: Properties of hadrons carrying open heavy flavor with charm or bottom quantum 
numbers C = +1 or B = +1, respectively 73 . Given are the valence quark content, the 


isospin, spin, and parity (I(J P )), the hadron mass, the most important decay modes together 
with their branching ratios (B.R.), and the decay length cr. Antihadrons with negative C or 
B quantum numbers are not listed explicitly. They decay into equivalent channels with charge 
conjugate decay products. 


Particle 

Quark 

content 

I(J P ) 

Mass (GeV/c 2 ) 

Decay mode 

B.R. (%) 

ct (/im) 

D+ 

cd 

l(o-) 

1.8696 ±0.0002 

K“7T + 7r + 
e + anything 
/j + anything 

9.13 ±0.19 

16.07 ±0.30 

17.6 ±3.2 

312 ±2 

D° 

cu 

1(0-) 

1.8648 ±0.0001 

K-7T+ 

e + anything 
ji + anything 

3.87 ±0.05 

6.49 ±0.11 

6.7 ±0.6 

123 ± 1 

D s + 

cs 

0(0-) 

1.9685 ±0.0003 

</>7T + 

e + anything 

4.5 ±0.4 

6.5 ±0.4 

150 ±2 

D *+ 

cd 

I(O-) 

2.0102 ±0.0001 

D° 7T+ 

67.7 ±0.5 

(2.1 ±0.5) x 10“ 6 

A c 

udc 

0(f) 

2.2865 ± 0.0001 

1 

+ 

5.0 ±1.3 

60 ±2 

B+ 

ub 

I(O-) 

5.2792 ± 0.0003 

J /ipK + 
l + vi anything 

0.1013 ±0.0034 
10.99 ± 0.28 

492 ±2 

B° 

db 

i(o-) 

5.2795 ± 0.0003 

iH’Ki 

l + vi anything 

0.0436 ± 0.0016 
10.33 ± 0.28 

455 ±2 

b hadrons 

sb 

0(0-) 

5.3663 ± 0.0006 

vn ' 

J/i/j anything 

0.14 ±0.05 

1.16 ±0.10 

441 ±8 


collisions from the pre-RHIC era are briefly summarized in Section [2} The main focus of this 
article is a review of the results on heavy-flavor yields, transverse momentum and rapidity 
distributions, as well as azimuthal distributions with respect to the reaction plane in heavy-ion 
collisions at RHIC (Section [3]) and at the LHC (Section [4]). Various theoretical approaches are 
confronted with the data and implications for the properties of the hot and dense QCD matter 
produced in such collisions are discussed. The article closes with a summary of the current 
status and an outlook to the future open heavy-flavor program at RHIC, LHC, and elsewhere 
in Section^ 


1.2 Open heavy flavor in pp and nucleus-nucleus collisions 


Open heavy-flavor production is experimentally accessible through the measurement of hadrons 
carrying charm or bottom quarks or through the measurement of their decay products. The 
most important properties of the relevant hadrons are summarized in Tab. [T 173 . 


In contrast to the light u, d, and s quarks, both the charm and the bottom quark are “heavy” 
with (bare) quark masses significantly exceeding the QCD scale parameter A qcd — 0.2 GeV 
(m c = 1.29loni GeV and m b = 4.19ig!o6 GeV 73 ). Thus, the production of a pair of a 


heavy quark and its antiquark in ultra-relativistic pp collisions proceeds exclusively through 
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initial hard partonic scattering processes and, therefore, can be treated theoretically within the 
framework of perturbative QCD. The unique feature of heavy-flavor production is that such a 
perturbative treatment is warranted for all momenta since the large quark mass introduces a 
hard scale even at zero momentum. This is in distinct contrast to gluon and light quark jets 
which can be treated perturbatively only at high p t . 

Consequently, the measurement of heavy-flavor production in pp collisions provides a cru¬ 
cial testing ground for pQCD. In the following, only hadro-production of heavy-flavor hadrons 
is addressed. A more detailed introduction, which also discusses photoproduction or the pro¬ 
duction of heavy quarks in e + e _ collisions, can be found elsewhere 74-76 . In this pQCD 


approach the basic production process of a hadron carrying heavy flavor can be factorized into 
three components: 


• the (non-perturbative) initial conditions. These are mainly determined by the fractional 
momenta, x, the interacting partons carry from the colliding hadrons. The x distributions 
of various partons inside hadrons have been studied mainly in deep inelastic scattering 
experiments. They depend on the squared energy-momentum transfer, Q 2 , between the 
two partons and they have been parametrized in the form of parton distribution functions 
(PDF). 

• the partonic (perturbative) scattering cross section. This can be calculated in perturbative 
QCD. At leading order (LO), the only processes contributing to heavy-flavor production 
are gluon fusion and quark anti-quark annihilation. At next-to-leading order (NLO) pro¬ 
cesses such as gluon splitting or flavor excitation have to be considered in addition. State 
of the art perturbative calculations of heavy-flavor production go even further in the sense 
that they calculate cross sections at fixed order with a next-to-leading-log resummation 
of higher orders in ai s . 

• the (non-perturbative) fragmentation of the heavy quarks into heavy-flavor hadrons. 
Here, two cases have to be distinguished: 


1 . formation of hadrons carrying open heavy flavor: in this case the heavy quark and 
antiquark individually fragment into hadrons. The energy of a heavy-flavor hadron 
with respect to the energy of the initially produced heavy quark is given by the 
fragmentation function, which is measured in e + e _ reactions and is assumed to be 
universal, i.e. the fragmentation of a heavy quark does not depend on the mechanism 
by which this quark was produced. 


2 . formation of quarkonium states: typically 1-2% of the produced heavy quark-anti¬ 
quark pairs form a bound quarkonium state instead of a pair of hadrons with open 
heavy flavor. While the physics of quarkonia is not discussed further in this article 
(see references 77-80 for reviews of this topic instead) it is important to mention 
that the measurement of the total charm and bottom yields would provide a natural 
reference for the investigation of charmonium and bottomonium yields. 


In the context of this review results from two modern NLO pQCD calculation approaches will 
be confronted with available open heavy-flavor hadro-production data: Fixed Order calculations 
with Next-to-Leading-Log resummation (FONLL) [81—831 and calculations within the General- 


Mass Variable-Flavor-Numbering Scheme (GM-VFNS) [84,85 













In addition to serving as a precision test of QCD, open heavy-flavor measurements in pp 
collisions provide a baseline for heavy-flavor studies in heavy-ion collisions, where the heavy 
quarks propagate through and interact with the produced hot and dense QCD matter. 

In the absence of nuclear effects, the heavy-flavor yields in nucleus-nucleus collisions would 
scale with the number of binary collisions as it is the case for all hard probes. Therefore, 
departures from binary scaling would indicate nuclear modifications of heavy-flavor observables. 
Two classes of such modifications have to be distinguished from each other, i.e. initial state 
effects and final state effects. 

Initial state effects result from the fact that the distributions of partons embedded in nuclei 
are different from the parton distribution functions in nucleons. Depending on x and Q 2 , 
different nuclear modifications of the parton distribution functions are observed. In the anti¬ 
shadowing region at medium x (around x ~ 0.1), the parton density in nuclei is larger than in 
nucleons. The shadowing region is characterized by a depletion of the parton density at low x 
(x < 10 -2 ) in nuclei relative to nucleons. Other effects, such as gluon saturation, may lead to 
reduced gluon densities at very low x as well. Reviews of nuclear parton distribution functions, 
on gluon saturation, and on related issues can be found elsewhere 


86 


Final state effects are in-medium modifications of heavy-flavor observables due to the pres¬ 
ence of hot and dense matter in ultra-relativistic heavy-ion collisions. The main focus of the 
open heavy-flavor programs at RHIC and at the LHC is the investigation of such modifications 
of open heavy-flavor observables, aiming to shed light on the properties of the hot QCD medium 
and the nature of parton-medium interactions. In turn, this requires a solid understanding of 
the initial state effects such that these and the final state effects can be disentangled. A number 
of final state effect has been predicted and is investigated experimentally. 

The total open heavy-flavor yield is sensitive to modifications of the initial parton densities 
but it is not expected to be changed significantly via final state effects. Because the charm 
and bottom quark masses are large enough that secondary, thermal heavy-flavor production 
mechanisms in the early, hot stages of a heavy-ion collision do not contribute significantly to 
the total heavy-flavor yield at RHIC. It is an interesting question whether thermal production 
of heavy quark-antiquark pairs can be observed in Pb+Pb collisions at the LHC but, most 
likely, this mechanism will not provide a large contribution to the total heavy-flavor yield even 
at LHC energies. 

A number of final state effects are expected to leave their footprint on the phase-space distri¬ 
butions of heavy-flavor hadrons. In particular, the interaction of partons with the medium 89 


is expected to be sensitive to the medium’s energy density via the mechanism of parton energy 
loss (see 90 for a recent review). This energy loss will reflect itself in a softening of the p t 


distributions of heavy-flavor hadrons or their decay products with respect to the same spectra 
in pp collisions. Prior to systematic heavy-flavor measurements at RHIC, it was generally ex¬ 
pected that the dominant energy loss mechanism of heavy quarks would be radiative energy loss 
via medium induced gluon radiation. Models incorporating this energy loss mechanism provide 
a reasonable description of the high-p t suppression of light-flavor hadrons at RHIC |91 -97 


However, in QCD the radiative energy loss is not the same for all partons. Gluons exhibit a 
larger color coupling factor than quarks such that the energy loss of quarks is expected to be 
smaller than that of gluons. Furthermore, with increasing quark mass induced gluon radiation 
at small angles with respect to the quark momentum vector is reduced owing to the so-called 
“dead cone effect” |98|. Radiative energy loss should therefore lead to a distinctive mass hier¬ 
archy of liigli-pt hadron suppression as quantified by the nuclear modification factor Raa- In 
the p t range up to about 10 GeV/c, where the heavy quark masses are sizable with respect to 
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their momenta, the energy loss is expected to decrease when going from light-flavor hadrons, 
which originate mostly from gluon or light-quark jets, to hadrons carrying charm and hadrons 
carrying bottom quarks: R AA < R AA < R\ A . 

It was pointed out, however, that collisional energy loss via elastic scattering should be 
comparable in magnitude to radiative energy loss over a wide kinematic region 99 . This 


additional contribution is more important for heavy than for light quarks and it is included in 
a number of model calculation which are confronted with data later in this review. 

Another class of models aiming at the description of heavy quarks in the QCD medium 
are Langevin-based transport models [ioo|[mi| in which a heavy quark is placed into a ther¬ 
mal medium. In such calculations the interaction is given by elastic collisions, in some ap¬ 
proaches [101] mediated by the excitation of D- or B-meson like resonant states in the medium. 
Heavy-quark diffusion and the damping of the initial non-equilibrium dynamics of heavy quarks 
in the medium are at the heart of these models, which gives indirect access to the viscosity to 
entropy density ratio of the QCD medium. 

heavy-flavor energy loss via repeated sequential formation 


In another approach 102,103 


of heavy-flavor hadrons and their dissociation in the medium is investigated. In this context, 
the interplay between the formation time of heavy-flavor hadrons and the expansion dynamics 
of hot QCD matter plays a key role. 

A completely different Ansatz |104| makes use of the AdS/CFT (anti de-Sitter-space/conformal- 
field theory) conjecture, which claims a correspondence between certain strongly coupled string 
theories and semiclassical gravitational physics. In this framework, heavy-quark energy loss is 
described in terms of string drag energy loss. 

In most of the approaches listed above it is assumed that the chemical composition of heavy- 
flavor hadrons is the same in pp and nucleus-nucleus collisions. However, in the light-quark 
sector an enhancement of baryons relative to mesons has been observed at intermediate p t 
in central heavy-ion collisions at RHIC 1105— 10T|. This becomes relevant as soon as heavy- 
flavor production is studied via the measurement of decay products of heavy-flavor hadrons, 
in particular electrons or muons from heavy-flavor semileptonic decays. For example, owing 
to a smaller semileptonic decay branching ratio and a softer decay lepton spectrum of A c 
baryon compared to D-meson decays, the Raa of electrons from charm hadrons decays would 
be smaller than one even without energy loss of charm in the medium 108,109 . Consequently, 


not only heavy-flavor mesons but also baryons have to be measured in exclusive decay channels 
in order to allow for unambiguous conclusions concerning in-medium modifications of heavy- 
flavor observables. Also for precision measurements of the total heavy-flavor production yields 
in pp and heavy-ion collisions it is mandatory to directly determine the contribution from 
heavy-flavor baryons. Unfortunately, until today heavy-flavor baryon yields and spectra could 
not be measured in heavy-ion collisions. 

Further insight into the interaction between heavy quarks and the hot QCD medium can be 
gained from the measurement of the azimuthal anisotropy of heavy-flavor hadron spectra with 
respect to the orientation of the reaction plane in heavy-ion collisions. In particular, the elliptic 
flow parameter v 2 is expected to be sensitive to the degree of thermalization of heavy quarks 
within the partonic medium. The measurement of a non-zero v 2 would certainly be indicative 
of a strong interaction of heavy quarks with the medium. In addition, heavy-flavor elliptic flow 
data are useful to benchmark models in which heavy-flavor hadrons are formed via the quark 
coalescence mechanism. Most of the theoretical approaches introduced above in the context of 
heavy-quark energy loss in the medium also provide predictions for v 2 such that a simultaneous 
measurement of v 2 and Raa provides a more stringent test of the underlying physics. 
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Figure 2: Equivalent pp production cross sections of neutral and charged D-mesons as functions 
of y/s in fixed-target experiments LO pQCD calculation with the PYTHIA event generator using 
different parton distribution functions are compared with the experimental data (left panel). 
Corresponding B-meson production cross sections from fixed target irp and pp collisions (right 
panel) (reprinted from Ref. l SI)- 


In summary, the value of heavy-flavor observables in hadronic collisions is twofold. First, in 
pp collisions perturbative QCD can be tested in a unique way and, thus, a solid experimental 
and theoretical reference can be prepared which is necessary for corresponding heavy-flavor 
measurements in nuclear collisions. Second, in collisions of the latter type the interaction 
of heavy charm and bottom quarks with hot and dense QCD matter is expected to be not 
only quantitatively but even qualitatively different from the interaction of light quark and 
gluon jets with this medium. Consequently, heavy-flavor measurements allow a complementary 
characterization of the parton-medium interactions and of the properties of strongly interacting 
matter, as it is produced in ultra-relativistic heavy-ion collisions, which can not be achieved 
otherwise. 


2 Open heavy-flavor hadro-production measurements in 
the pre-RHIC era 

2.1 Open heavy flavor from fixed target experiments 

The discovery of hadrons carrying open charm or bottom was preceded in both cases by the 
discovery of charmonium or bottomonium, respectively. The J/^ charmonium state was discov¬ 
ered in the so-called “November revolution” simultaneously in p-Be collisions at the Brookhaven 
AGS (l 101 and in e + e _ collisions at Stanford’s storage ring SPEAR 111 . The T bottomonium 


state was discovered in p-nucleus collisions at FNAL (112]. Experimentally, the advantage of the 
J /xp and the T compared to open heavy-flavor hadrons is related to the fact that the quantum 
numbers of these quarkonia coincide with those of the photon. Therefore, they can be produced 
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resonantly in e + e _ colliders for precision measurements of their properties. More important in 
the context of heavy-ion collisions, the J/r/ 1 and T can decay into dilepton pairs giving rise to 
the unique experimental signature of a high p t dilepton pair, which can be triggered on. Both 
neutral |113| and charged 114| open-charm D mesons were first observed at the SLAC SPEAR 
storage ring in purely hadronic decays. First evidence for open-bottom production was found 
with the CLEO 1115 and CUSB 1116; detectors at the Cornell Electron Storage Ring (CESR) 
via the measurement of electrons from semileptonic B-meson decays. 

In contrast to heavy quarkonia, in particular the J /ip meson, which played a crucial rule in 
the fixed-target heavy-ion program at the CERN SPS, it was only with the advent of high energy 
nucleus-nucleus colliders, i.e. RHIC and LHC, that open heavy-flavor hadrons could take their 
role as unique probes for hot QCD matter. For many years open charm and bottom production 
was in the focus of elementary particle physics only, studied mostly via photo-production or in 
e + e _ experiments. 

Although single electrons from heavy-flavor hadron decays were first observed in pp collisions 
at the CERN ISR collider at y/s = 52.7 GeV 117 , actually before the discovery of charm, 


subsequent studies of open heavy flavor hadro-production were conducted mainly in fixed-target 
experiments at y/s < 50 GeV employing pion and proton beams at the CERN SPS, at FNAL, 
and at DESY. Most of these measurements, which are reviewed in detail elsewhere 118 , were 


specifically designed for heavy-flavor measurements. As a key feature they included detectors 
with high spatial resolution in the target region in order to separate the primary collision 
vertex from secondary heavy-flavor decay vertices. Bubble chambers, emulsions, and silicon 
tracking telescopes were used for this purpose. Large statistical samples have been collected 
with pion beams only. The proton beam experiments accumulated much less data, with the 
~ 300 neutral and charged D-meson events collected by E769 representing the largest statistics 
data sample 11191. It is interesting to note that not only E769 but also most of the other 


experiments used nuclear targets, notably Be, Al, Cu, and W. To derive equivalent production 
cross sections for pp (or 7rp) reactions a linear dependence of the cross section on the mass 
number A of the target nucleus was assumed, i.e. cr p A, 7 rA = A • <7 PP;7rp . In Fig. [2] the derived D- 
and B-meson pair production cross sections in pp or irp collisions, extrapolated to 47t from the 
phase space covered in the various experiments, are shown as a function of y/s. The leading- 
order pQCD inspired event generator PYTHIA 11201 can be used to describe describe heavy- 
flavor production in the fixed-target energy regime as was demonstrated in Ref. 11211. More 


detailed and systematic studies have substantiated these findings as summarized in Ref. [! 18 


Calculations with PYTHIA employing different parton distribution functions are compared 
with charged and neutral D-meson production cross sections from pp collisions in the left panel 
of Fig. [2| The fact that higher-order diagrams are absent in PYTHIA is compensated by 
scaling the calculated cross sections up by an empirical factor, the so-called K-factor. This 
approach is justified only under the assumption that the kinematical distributions of heavy 
quarks originating from leading order and higher order processes are similar. Depending on the 
choice of other model parameters in PYTHIA, K-factors in the range 2.5 - 4.5 are necessary 
to obtain a satisfactory agreement of the calculation with charm production data. The bottom 
production cross section is smaller than the charm production cross section by almost three 
orders of magnitude in this y/s range as shown in the right panel of Fig. [2j The experimental 
uncertainties are large and some of the measurements, e.g. the E771 and E789 results obtained 
at the same energy, are inconsistent with each other. It is, however, remarkable that PYTHIA 
calculations, using a K-factor of 2 for bottom production, do not only describe the low energy 
bottom production cross sections shown in Fig. [2] within a factor of 2, but they are also in 
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Figure 3: Charm hadro-production cross section in the forward hemisphere (xp > 0) as a 
function of y/s in comparison with the \fs dependence of a corresponding PYTHIA calculation. 
The absolute normalization of the PYTHIA curve was fitted to the experimental data (reprinted 
from Ref. 1124 with kind permission from Springer Science and Business Media). 


reasonable agreement with bottom production data measured in pp collisions at the Tevatron 
Collider at y/s close to 2 TeV as discussed in Ref. 118 . The same systematic analysis of 


heavy-flavor hadro-production data has demonstrated that for y/s values accessible at fixed- 
target machines not only the open heavy-flavor absolute production cross sections but also the 
D-meson kinematical distributions and DD pair correlations can be described reasonably well 
with PYTHIA. 


While in the heavy-ion program at the CERN SPS open heavy-flavor production was never 
measured directly, the contribution of simultaneous semimuonic decays of correlated DD pairs to 
the dimuon continuum was carefully studied. Particular attention was given to the intermediate 
mass region (IMR) of the dimuon mass distribution, i.e. the region between the <fi and the 
J/ip. In addition to open charm decays, the only relevant sources for dimuons in the IMR are 
the Drell-Yan process and, potentially, thermal dimuon production in the hot QCD matter in 
nucleus-nucleus collisions. Intermediate mass dimuon production was first measured at the SPS 
with spectrometers that did not allow for a separation of prompt dimuons from the primary 
collision vertex and dimuons from displaced, secondary decay vertices characteristic for open 
charm decays. The NA38 1221 and HELIOS-3 Collaborations 123 measured dimuons in p- 
W, S-U, and S-W collisions at a beam energy of 200 GeV per nucleon. Later, NA38/NA50 
investigated p-A collisions, using Al, Cu, Ag, and W as nuclear targets at a beam energy of 
450 GeV, and Pb+Pb collisions at 158 GeV per nucleon [124 . 


The correlated, opposite-sign IMR dimuon spectra measured by the three experiments in 
p-A collisions are in good agreement with the sum of the expected Drell-Yan and open charm 
contributions. In fact, NA50 determined the ratio of dimuons from open charm decays to 
Drell-Yan dimuons via fits to the mass spectra, and used PYTHIA to infer indirectly the total 
charm production cross section in pp collisions, assuming that in p-A collisions the total open 
charm cross section scales with the nuclear target mass number. The resulting charm cross 
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Figure 4: Dimuon yield measured in peripheral ((N paxt ) = 110, left panel) and central ((lV part ) = 
381, middle panel) Pb+Pb collisions at a beam energy of 158 GeV per nucleon in comparison 
with the contributions from expected sources. The sum (solid) of Drell Yan dimuons (dashed), 
open charm (dotted), and the dimuon decay of charmonia (dashed-dotted) is not sufficient 
to saturate the measured IMR dimuon yields. Excess relative to the expected open charm 
contribution as a function of the number of participants for various collision systems (right 
panel) (reprinted from Ref. 1124 with kind permission from Springer Science and Business 
Media). 


section is shown in Fig. [3] in comparison with direct measurements of the charm production 
cross section and as a function of yfs together with a PYTHIA calculation similar to the one 
discussed above. The cross section obtained indirectly from the IMR dimuon spectra in p-A 
collisions agrees reasonably well with the yfs dependence obtained from direct measurements. 

In contrast, in heavy-ion collisions the measured yields of IMR dimuons exceeded the ex¬ 
pectation from Drell-Yan and open charm decay contributions as demonstrated in the two left 
panels of Fig. |4} The excess, which is clearly visible in the IMR, is more pronounced in central 
than in peripheral Pb+Pb collisions. The centrality dependence is quantified in the right panel 
of Fig. [4] which shows the excess relative to the expected dimuon yield from open charm decays 
as a function of the number of participants as measured by NA38 and NA50 in various collision 
systems. 


Among several other potential explanations of the observed IMR dimuon excess two inter¬ 
pretations were discussed as the most likely scenarios, i.e. thermal production as an additional 
source of IMR dimuons or an increase of the open charm production cross section per nucleon in 
heavy-ion collisions compared to p-nucleus collisions. The clarification of the origin of the IMR 
dimuon excess was one of the main goals of the NA60 experiment. The NA60 setup combined 
the muon spectrometer (MS) previously used by the NA38 and NA50 experiments 125 with a 


radiation hard silicon vertex tracker (VT) located in a dipole magnet providing a magnetic field 
of up to 2.5 T 1261. The VT was positioned in between the interaction targets and the upstream 
hadron absorber in front of the MS. Compared to the previous experiments, the unique feature 
of this setup was the capability to match tracks of charged particles reconstructed in the VT 
and the MS both in position AND in momentum space, thus combining the high-rate capability 
and selective dimuon trigger of the MS with the excellent position and improved momentum 
resolution provided by the VT. This made it possible to investigate dimuon production in In-In 
collisions at a beam energy of 158 GeV per nucleon with an unprecedented precision. The 
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presence of an IMR dimuon excess consistent with the earlier observations in Pb+Pb collisions 
was confirmed with the NA60 experiment in the In-In collision system 1127,128). Furthermore, 
it was proven that the excess yield does not originate from a displaced, secondary vertex but 
is consistent with prompt dimuon emission from the primary collision vertex [127[|128| , thus 
ruling out an enhanced charm production cross section as source responsible for the excess. In¬ 
stead, the measured properties of the excess dimuons suggest that they are related to thermal 
radiation from the hot QCD medium. 

In summary, open heavy-flavor hadro-production in fixed target experiments was observed 
to be consistent with expectations from perturbative QCD. This includes data from fixed target 
heavy-ion collisions, which to date are available only from the SPS. Since open heavy-flavor 
hadrons were never directly measured in fixed-target heavy-ion collisions it is not precluded 
that such measurements can provide useful and unique information about the properties of 
strongly interacting matter in the future. In fact, precision measurements of charm production 
in nucleus-nucleus collisions close to the production threshold is one of the key goals of the 


Compressed Baryonic Matter (CBM) experiment 129 at the future FAIR facility. 


2.2 Open heavy flavor from proton/antiproton colliders 

Before the startup of RHIC, open heavy-flavor production was investigated in elementary col¬ 
lisions at other hadron colliders giving access to much higher y/s compared to previous fixed 
target experiments. For the discussion of medium modifications of open heavy-flavor observ¬ 
ables in heavy-ion collisions it is useful to investigate whether pQCD calculations provide a 
theoretical baseline for heavy flavor hadro-production at high collision energies. 

The very first open heavy flavor hadro-production measurements at a collider date back to 
the experimental program at the CERN Intersecting Storage Rings (ISR) which is reviewed 
in detail elsewhere (l30l - Electrons from semileptonic charm hadron decays were observed 
for the first time in the early 1970s in pp collisions at y/s = 52.7 GeV at the ISR 1117 


However, at the time of these experiments the origin of single electrons measured in the range 
1.6 < p t < 4.7 GeV/c was not known as the charm quark was not discovered yet. Further 
measurements of single electrons and lepton pairs (ee and ep) as well as full reconstructions of 
hadronic decays of open charm hadrons indicated substantial charm production cross sections 
at ISR energies. 

Bottom production at hadron colliders was addressed first in the late 1980s at the CERN 
SppS. The UA1 Collaboration investigated bottom production at high p t (p t > mb, where 
mb is the bottom quark mass) in pp collisions at y/s = 546 GeV and 630 GeV [131,132]. 
For this purpose, events containing single high-p t muons (10 < Pt <40 GeV/c) together 
with a jet or events containing dimuon pairs (2 < m^ < 35 GeV/c 2 ) were selected. Applying 
isolation cuts on the single muons and dimuon pairs, respectively, the contributions from bottom 
hadron decays to the data could be determined statistically. The resulting p t -differential bottom 
production cross sections from UAl were compared with next-to-leading-order, i.e. order a^, 


pQCD calculations which were state of the art calculations at that time 133-136 . In general, 


a good agreement between data and pQCD calculations was observed within the substantial 
experimental and theoretical uncertainties. The measurement of open-charm production was 
difficult at the SppS. At low p t , no selective triggers were available, and at high p t it was 
difficult to separate charm from bottom production. The D* + meson was the only charmed 
hadron which could be directly reconstructed with the UAl experiment in hadronic decays at 


high p t in jets 137 . With the UA2 experiment single electrons from semileptonic open charm 
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Figure 5: /^-differential cross section of J/b mesons originating from the decay of bottom 
hadrons integrated over \y\ < 0.6 as measured with the CDF experiment in pp collisions at yfs 
= 1.96 TeV (left panel) and the corresponding bottom hadron p t differential cross section de¬ 
termined from the J/'i/j data (right panel). The plotted J/A cross section includes the branching 
ratio of Jmesons in the dimuon decay channel. FONLL pQCD calculations are compared 
with the data. The solid lines show the central theoretical values and the dashed lines indicate 
the theoretical uncertainties (reprinted with permission from Ref. [150 
the American Physical Society). 
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hadron decays were measured in the transverse momentum range 0.5 < p t < 2.0 GeV/c fl38 


Only ten such events were observed in this experiment which was the first to use a RICH 
detector for electron identification at a collider. A detailed review of heavy-quark production 
at the CERN SppS can be found elsewhere 1139 


Open heavy-flavor production measurements, mainly focusing again on bottom production, 
were conducted in the 1990s at the Fermilab Tevatron with the CDF and DO experiments in 
pp collisions first at yfs = 1.8 TeV (Run-I) and later at yfs = 1.96 TeV (Run-II). The CDF 
Collaboration reported the first full reconstruction of B mesons (using the decay B ± —>■ J/AK ± ) 
at high p t (above 9 GeV/c) and at midrapidity (|y| < 1.0) at a hadron collider 1401. Further 
bottom measurements with the CDF experiment employing electrons El and muons [142] 
from semileptonic bottom decays as well as B° mesons in the exclusive decay channel B° 
J/V>K*(892) [M3; seemed to have established that bottom production cross sections measured in 
the 1.8 TeV Run-I at the Tevatron exceeded the expectation from pQCD calculations. While the 
first bottom-production measurements with the DO experiment, based on single high-p t muons 
at midrapidity, were in reasonable agreement with NLO pQCD calculations 144 , later higher 


statistics measurements of single muons and dimuons j 145,146 , and J/p mesons from displaced, 
secondary vertices at midrapidity 147 followed the trend suggested by the CDF experiment 
that bottom was produced in excess of the yields expected from NLO pQCD calculations 
of heavy quark production. This apparent discrepancy between measured bottom production 
cross sections in pp collisions and corresponding NLO pQCD calculations posed a long standing 
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Figure 6: /^-differential cross sections of various D-meson species measured with the CDF ex¬ 
periment in pp collisions at y/s = 1.96 TeV in the rapidity range \y\ < 1. The solid curves are 
theoretical predictions from FONLL pQCD calculations with the shaded bands indicating the 
theoretical uncertainties. Predictions for the D+ cross section are not available from FONLL. 
Ratios of the measured cross sections and theoretical predictions from GM-VFNS pQCD cal¬ 
culations as functions of p t (reprinted with permission from Ref. 1511 (upper four panels) and 
Ref. [85] (lower four panels); Copyrights (2003, 2006) by the American Physical Society). 


puzzle. 

As it is documented in detail elsewhere 148,149 the Tevatron bottom puzzle was finally 
resolved with a more careful comparison of experimental and theoretical observables and with 
the availability of much improved data from the Tevatron Run-II at y/s = 1.96 TeV. In this 
context, the most important upgrades of the CDF and DO experiments were related to the trig¬ 
ger systems selecting the events to be recorded for further analysis. Both experiments included 
silicon vertex detectors in their trigger schemes in Run-II, which allowed them to trigger on the 
displaced secondary vertices of open heavy-flavor hadron decays. Due to bandwidth limitations 
the DO displaced vertex trigger could operate only on preselected events that contained an ad¬ 
ditional muon. The high-bandwidth displaced vertex trigger of the CDF experiment provided 
unique access to the measurement of fully hadronic open heavy-flavor hadron decays. Employ¬ 
ing this trigger, with the CDF experiment the /^-differential production cross section of J /ip 
from open bottom hadron decays as well as the underlying bottom hadron production cross 
section was measured down to low p t in pp collisions at y/s = 1.96 TeV 150 . As shown in 


Fig.@ this measurement is in excellent agreement with FONLL pQCD calculations as are all 
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further open bottom hadron production results from Run-II at the Tevatron. 


The CDF silicon vertex trigger did not only allow for a measurement of open bottom produc¬ 
tion in a unique way at the Tevatron Run-II but it gave also access to open charm production. 
The /^-differential production cross sections of open charm mesons were measured 151 , at 


high p t and at midrapidity ( \y\ < 1), in the fully reconstructed hadronic decays D° —> K~ 7r + , 
D*+ —)-D° 7r + , D + —t RT _ 7r + 7r + , D+ —)■ <fi 7r + , and their charge conjugates. The measured cross 
sections are shown in the upper four panels of Fig. [6| Results from FONLL pQCD calcula¬ 
tions 181,821 are plotted in comparison with the data except for the case of the D+ which is 
not considered explicitly in FONLL. While the calculations agree with the data within the 
experimental and theoretical uncertainties, a general trend is observed that the measured cross 
sections he at the upper edge of the calculated cross sections. This is further demonstrated 
in the lower four panels of Fig. [6] which show the ratios of the measured cross sections and 
corresponding pQCD calculations employing the GM-VFNS approach |85], 


In summary, open heavy-flavor hadro-production measurements conducted at pp and pp 
colliders before the startup of RHIC can be reasonably well described by state of the art pQCD 
calculation. Given the significant theoretical uncertainties and the fact that open heavy-flavor 
measurements have been proven to be difficult, it is crucial to establish an experimental open 
heavy-flavor production baseline from elementary pp collisions at RHIC, at the LHC, and at 
future accelerators investigating heavy-flavor hadro-production. Furthermore, such experimen¬ 
tal baselines have to be confronted with pQCD calculations before open heavy-flavor hadrons 
can be used with confidence as probes for the medium produced in ultra-relativistic heavy-ion 
collisions. 


3 Open heavy-flavor measurements at RHIC 

3.1 Experiments studying heavy-flavor production at RHIC 

RHIC delivers heavy-ion beams since the year 2000. Two of the four heavy-ion experiments at 
RHIC, PHENIX and STAR, have conducted systematic open heavy-flavor measurements since 
then, mainly in pp and Au+Au collisions at ^/snn = 200 GeV, where substantial data sets 
have been recorded. In addition, data from d+Au collisions, which give access to potential 
cold nuclear matter effects, and from Cu+Cu collisions, which help in the investigation of the 
system size dependencies of open heavy-flavor observables, were investigated. 

The results reviewed in this article originate from published measurements with the PHENIX 
and STAR experiments in their configurations prior to the years 2008 and 2010, respectively. 
Since then, both PHENIX and STAR have embarked on ambitious detector upgrade programs. 
Key components in both upgrade projects are precision vertex spectrometers which will allow 
the separation of tracks originating from displaced, secondary open heavy-flavor hadron decay 
vertices from primary tracks pointing back to the collision vertex. These improvements of the 
experimental setups go hand in hand with an upgrade of the RHIC accelerator facility aiming 
at an increase of the original Au+Au luminosity by a factor of five. The open heavy-flavor 
physics potential of the detector and accelerator upgrades is discussed in detail elsewhere 1152 . 
In the following, the experimental status quo of the PHENIX and STAR experiments as of the 
years 2007 and 2009, respectively, is briefly summarized. 
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Figure 7: Schematic beam view (left panel) and side view (right panel) of the PHENIX appa¬ 
ratus in the year 2005. Detectors relevant for the heavy-flavor measurements discussed here are 
the drift chambers (DC) and the pad chambers (PC) for tracking, and the Cherenkov detectors 
(RICH) and the electromagnetic calorimeters (PbGl, PbSc) for electron identification in the 
central arms. The relevant detectors in the forward and backward muon spectrometers are the 
muon trackers (MuTr) and the muon identifiers (MuID). Global beam-beam counters (BBC) 
and zero-degree calorimeters (ZDC) are used for event characterization (from Ref. 154]). 


3.1.1 The PHENIX approach to open heavy-flavor measurements 

The PHENIX experiment comprises four spectrometers and a set of global detectors for event 
characterization. The two central spectrometer arms, each covering 90° in azimuth, allow the 
measurement of electrons, photons, and charged hadrons in the pseudorapidity range \rj\ < 0.35. 
The momenta of charged particles are measured in the central arms using drift chambers and 
pad chambers located in a magnetic field which is oriented in parallel to the beam axis around 
the interaction vertex. For particle identification, time-of-flight detectors, RICH detectors, and 
electromagnetic calorimeters are employed in the central spectrometers. In the forward and 
backward regions, two muon spectrometers cover the pseudorapidity ranges —2.2 < q < —1.2 
and 1.2 < q < 2.4 over the full azimuth. Both muon spectrometers comprise a muon tracker 
with three stations of multi-plane drift chambers placed in a radial magnetic field, and a 
muon identifier consisting of alternating layers of streamer tubes and steel absorbers. As 
global detectors, beam-beam counters and zero-degree calorimeters are located at forward and 
backward rapidity, and they serve a variety of purposes. With these detectors the collision 
time and vertex position along the beam direction is determined. Furthermore, the first level, 
minimum bias trigger signals for collisions are derived from the BBC and the ZDC. In addition, 
they are used to determine the collision centrality and the azimuthal orientation of the reaction 
place in nucleus-nucleus collisions. Schematic beam and side views of the PHENIX experimental 
setup are shown in Fig. [7| A detailed description of the PHENIX apparatus can be found 
elsewhere 


153 


In the PHENIX configuration discussed here open heavy-flavor hadrons are not recon¬ 
structed exclusively via their hadronic decays. Instead, electrons from heavy-flavor hadron 
decays are measured with the central spectrometer arms around mid-rapidity in pp and nucleus- 
nucleus collisions. Both the transverse momentum spectra of single electrons and the dielectron 
mass spectra are inspected for contributions from semileptonic heavy-flavor hadron decays. Cor- 
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respondingly, with the forward/backward muon spectrometers single muons from heavy-flavor 
hadron decays are measured in pp and nucleus-nucleus collisions. Both in the electron and in 
the muon measurement, the subtraction of background from sources other than heavy-flavor 
hadron decays poses the most important experimental challenge. 


In case of the electron measurements explained in detail in Ref. 155 , remaining background 
from misidentified hadrons does not play a significant role since the combination of Cherenkov 
detectors and electromagnetic calorimetry guarantees the selection of a rather clean electron 
candidate sample. From the measured inclusive single electron p t distributions electron back¬ 
ground from non-heavy-flavor sources has to be subtracted. The most important components 
of this so-called background cocktail, which is calculated using a MonteCarlo hadron decay 
generator, are Dalitz decays of the 7 r° meson and the conversion of photons in detector mate¬ 
rial, where the majority of photons originates from the 7 r° —y 77 decay. It is important to note 
that, in the PHENIX experiment, special emphasis is put on the minimization of the material 
budget in the detector aperture, thus minimizing the contribution from photon conversions to 
the measured electron spectra. With a relevant material budget of less than 1% of a radiation 
length, yield ratios of electrons from photon conversion to electrons from light meson Dalitz 
decays significantly below one are achieved in the PHENIX experiment. Similar Dalitz decay 
and conversion contributions from other light mesons, notably the rj meson, are also considered 
in the electron background cocktail. Furthermore, dielectron decays of light vector mesons (p, 
u, 0 ) are included although these contributions played a negligible role only. At low p t , a 
minor contribution from weak kaon decays (K e3 ) has to be added to the background cocktail. 
Towards high electron p t (above 5 GeV/c) hard quark-gluon Compton scattering processes, 
contributing to the electron background via virtual photon production and via the conversion 
of direct real photons in material, play an increasingly important role. In this high p t region 
decays of heavy quarkonia (J /ip, T), and electrons from the Drell-Yan process can not be ig¬ 
nored either, in fact these sources represent about one third of the background electron yield 
in the range 5 < p t < 10 GeV/c in pp and Au+Au collisions at ^/snn = 200 GeV at RHIC. 
In the first electron measurements with PHENIX, which did not reach high electron p t with 
good statistics, contributions from J/r/> and T decays as well as from the Drell-Yan process 
were not included in the electron background cocktail 157-162 . Only with the availability 


of more precise data these sources were added to the cocktail 1155j. The background cocktail 


from these sources was subtracted from the inclusive electron sample and the difference was 
identified with electrons from open heavy-flavor hadron decays. A basic requirement for this 
cocktail subtraction method was the availability of the yields and phase space distributions of 
all relevant background electron sources, ideally from a measurement in the same apparatus. 
In addition, the applicability of the cocktail background subtraction method was in general 
limited to the electron p t range in which the ratio of signal electrons from heavy-flavor hadron 
decays to background electrons from other sources was larger than the combined systematic un¬ 
certainty of the inclusive electron measurement and the calculated background cocktail. Since 
the signal-to-background ratio decreases with decreasing electron p t , at low p t (typically below 
2 GeV/c in the PHENIX case) it became advantageous to switch from the cocktail subtraction 
to the so-called converter subtraction method to determine the yield of electrons from heavy- 
flavor hadron decays. In this complementary approach, additional converter material with a 
well defined thickness was added to the experimental setup for a fraction of the running time of 
the experiment. The converter multiplied the yield of “photonic” electrons, i.e. electrons from 
photon conversions and from Dalitz decays of light mesons, by a known factor, thus providing 
a measurement of this background contribution which was dominant at low electron p t . A 
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detailed description of the measurement of single electrons from heavy-flavor hadron decays 
with the PHENIX experiment using the cocktail and converter subtraction techniques can be 
found elsewhere |l55| . 

A disadvantage of single lepton measurements is the fact that contributions from charm and 
bottom hadron decays can not be separated from each other without additional experimental 
information. Electron-hadron correlation measurements can provide this missing information 
because electrons from heavy-flavor hadron decays are always emitted together with hadrons. 
Since the correlations are not the same in charm and bottom hadron decays they can be used to 
disentangle the corresponding relative contributions to the electron spectra from heavy-flavor 
hadron decays |I63| . 

The contributions to the electron-positron mass spectrum from simultaneous semielectronic 
decays of correlated pairs of hadrons carrying heavy quarks and antiquarks, respectively, was 
determined with the PHENIX experiment in a cocktail approach as well 1164,165 . Except 
for photon conversions, the same background sources as listed above for the case of single 
electrons were considered for the calculation of a background cocktail. After the subtraction 
of this cocktail background contributions from correlated heavy-flavor hadrons decays became 
accessible. 

Also for the measurement of muons from semileptonic heavy-flavor hadron decays in the 
forward and backward muon spectrometers of the PHENIX experiment the subtraction of 
background was crucial 1156,166 


In this case, three background sources were identified as 
the most important ones. First, decay muons from light hadrons decaying before reaching the 
first absorber material could not be distinguished experimentally from muons from heavy-flavor 
hadron decays. Second, punch-through hadrons produced at the collision vertex could penetrate 
through the MuID layers and, therefore, were included in the inclusive muon candidate sample. 
Third, muons from the decay of hadrons, which penetrated the muon arm front hadron absorber 
and decayed in the muon tracker volume, contributed to the inclusive muon candidate sample. 
The overall background from light hadron sources was quantified via MonteCarlo simulations. 
The contribution from decay muons could be subtracted statistically. This was possible because 
the yield of light hadron decay muons exhibited a pronounced vertex position dependence along 
the beam direction which was not the case for the other sources. Again, a detailed description 
of the experimental techniques used to isolate the contribution from open heavy-flavor hadron 


decays to the measured single muon spectra can be found elsewhere 156,166 


3.1.2 The STAR approach to open heavy-flavor measurements 

The subdetectors of the STAR experiment relevant in this context were housed inside a solenoid 
magnet providing a field of 0.5 T. Charged particle tracking relied mostly on a large azimuthally 
symmetric time projection chamber (TPC) which covered the pseudorapidity range \r]\ < 1.8. 
For tracking close to the interaction point and for vertexing the TPC originally was augmented 
by a silicon vertex tracker (SVT). This device, however, could not be used for a clean separation 
of decay vertices from open heavy-flavor hadron decays from the primary collision vertex and, 
hence, it was removed from the STAR setup prior to the RHIC run in the year 2008 to reduce 
the material budget in front of the TPC inner field cage by a factor of ~ 10. For particle 
identification, the specific energy loss measured in the TPC could be combined with velocity 
information from a Time-of-Flight barrel which employed multigap resistive plate technology. 
For electron identification, a barrel electromagnetic calorimeter (EMC) was available in addi¬ 
tion. This device included a preshower and a shower maximum detector. As in the PHENIX 
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Figure 8: Schematic view of the STAR detector. Charged particles tracks are reconstructed 
around mid-rapidity in a large volume time projection chamber (TPC), which also contributes 
to the particle identification via the measurement of the charged particles’ specific energy loss 
dE/dx. For additional particle identification the TPC is surrounded by a Time-of-Flight barrel 
(ToF) and an electromagnetic calorimeter (EMC) (from Ref. 1168]). 


case, beam-beam counters and zero degree calorimeters could be used to provide global event 
characteristics. A schematic view of the STAR experimental setup is shown in Fig. [8] and a 
detailed description of the STAR apparatus can be found elsewhere |l67| . 


It is only with the STAR experiment that heavy-flavor hadrons have been fully recon¬ 
structed at RHIC via an invariant mass analysis of their hadronic decay products. The main 
difficulty of such measurements is the substantial combinatorial background from other uncorre¬ 
lated hadrons. With increasing charged particle multiplicity the resulting signal-to-background 
ratio becomes worse such that the difficulty of such a measurement increases from pp to 
proton/deuteron-nucleus to nucleus-nucleus collisions in which signal-to-background ratios of 
typically 1:1000 had to be coped with in the STAR experiment. With STAR, D° mesons have 
been reconstructed in their K-7r decay channel in d+Au 11691 and, recently, in pp collisions 1170 


A corresponding signal was observed in Au+Au collisions at low p t as well 171 . In addition 


to the D° also the D* + meson has been reconstructed via its D°7 t + decay in pp 1172 but not in 
nuclear collisions. The only possibility to increase the signal-to-background ratio is the clean 
separation of the secondary decay vertex from the primary collision vertex which has not been 
achieved with the experimental setup of the STAR experiment discussed here. 

Heavy-flavor measurements in Au+Au collisions were conducted with the STAR experiment 
mainly via the semielectronic decays of heavy-flavor hadrons |l73 , which were also investigated 
for reference in the pp collision system 173,174 . Electrons were identified via their specific 


energy loss dE/dx in the STAR TPC, their velocity as measured with the Time-of-Flight detec¬ 
tor system, and their response in the electromagnetic calorimeter. A remaining contamination 
from misidentified hadrons in the electron sample was quantified using the TPC dE/dx mea- 
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surement, and it was subtracted statistically. As in the PHENIX case, the spectra of electrons 
from open heavy-flavor decays were obtained via the subtraction of electron background from 
the measured inclusive electron spectra. Originally, the electron background in the STAR ex¬ 
periment was completely dominated by photon conversions in detector material, in particular 
in the SVT. With the removal of this detector for the 2008 run at RHIC, the relevant material 
budget was reduced by a factor sMO to just above 1% of a radiation length, which was close 
to the material budget in the PHENIX setup. The background from photon conversions and 
Dalitz decays was, however, determined distinctively different within STAR compared to the 
PHENIX cocktail and converter subtraction approaches. Electron-positron pairs from the rel¬ 
evant background sources have a very small invariant mass, below the 7 r° mass. In contrast to 
the PHENIX experiment, the large acceptance STAR apparatus provided a reasonably large 
efficiency (~ 60%) for the reconstruction of these low mass pairs, such that this contribution to 
the background could be measured and subtracted from the inclusive electron spectrum. Not all 
of this so-called “photonic background” could be identified with this technique because one of 
the electrons might have been outside of the detector acceptance or had a momentum to small 
to be measured. In this case both electrons from the pair were not reconstructed as part of the 
photonic background. This inefficiency had to be determined precisely through simulations. 
Contributions from other background sources, e.g. dielectron decays of vector mesons which 
were not measured directly, were calculated in a cocktail approach. The difference between 
the measured inclusive electron spectrum and the resulting total background from photonic 
and other sources was identified with the spectrum of electrons from open heavy-flavor hadron 
decays. 

In pp collisions the azimuthal angular correlations of electrons from open heavy-flavor 
hadron decays with hadrons or D° mesons was investigated with the STAR experiment 
These procedures allowed a statistical separation of electrons from charm and bottom hadron 
decays, similar to the approach chosen by PHENIX. 

While the STAR experiment does not include a dedicated muon spectrometer muons could 
be identified at very low p t (close to 0.2 GeV/c) via their time of flight and their ionization 
energy loss in the TPC HZil- Muons from heavy-flavor hadron decays could be separated sta¬ 
tistically from the dominant background from weak pion and kaon decays via the measurement 
of the distance of closest approach of the tracks to the primary collision vertex. The resulting 
yield of low-p t muons from heavy-flavor hadron decays provided some sensitivity to the total 
open charm cross section [171]. 

3.2 Reference measurements in pp collisions 

3.2.1 Overview 

The open heavy-flavor measurements with the PHENIX and STAR experiments in pp collisions 
at y/s = 200 GeV at RHIC aim at two goals. First, these data provide a crucial experimental 
reference for corresponding measurements in nucleus-nucleus collisions. Such a reference is a 
key ingredient for the investigation and interpretation of effects due to the presence of a hot and 
dense medium produced in high-energy nucleus-nucleus collisions. Second, these data provide 
an important testing ground for perturbative QCD calculations which, due to the large masses 
of the charm and bottom quarks, should be able to predict open heavy-flavor observables even 
at low p t - 

As will be demonstrated in this section, predictions from NLO pQCD calculations within 
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Figure 9: Raw D* + candidate signal from minimum bias pp collisions at ^/s = 200 GeV as 
measured with the STAR experiment, (left panel) Results from two methods to quantify the 
combinatorial background are shown in addition (see text for details), cc production cross 
section as deduced from the measured D° and D* + cross section (right panel). A power-law 
fit and the p t differential charm cross section from a FONLL pQCD calculation are shown 
for comparison (reprinted with permission from Ref. [170] ; Copyright (2012) by the American 
Physical Society). 


the FONLL approach |81,82 are in reasonable agreement with essentially all observables related 


to the production of open heavy-flavor hadrons in pp collisions at RHIC. 


3.2.2 Hadronic heavy-flavor hadron decays 

The only full reconstruction of hadronic decays of heavy-flavor hadrons in pp collisions at 
200 GeV at RHIC was the measurement of the D° — > K _ 7t + decay in the range 0.6 < p t < 
2.0 GeV/c and the D* + — > D°7t + decay in the range 2.0 < p t < 6.0 GeV/c with the STAR 
experiment [170]. As an example, the left panel of Fig. [9] indicates the quality of the D* + 
reconstruction. For minimum bias events, the distribution of the invariant mass difference 
AM = M(K7T7 t) — M(Ktt) is shown for “right sign” combinations (K =F 7r ± 7r ± ), which contain the 
D* + signal, for “wrong sign” combinations (K ± 7r =F 7r ± ), which are combinatorial background, 
and for D° “side band” combinations, which select K7 t pairs that are not consistent with the D° 
meson mass. Both the “wrong sign” and the “side band” method reproduce the combinatorial 
background well. After corrections the measured p t differential production cross sections for 
D° and D* + mesons in pp collisions at y/s = 200 GeV were divided by the charm quark 
fragmentation ratios 0.565 ± 0.032 for c —> D° and 0.224 ± 0.028 for c —> D* + , respectively, 
to convert the meson cross sections to the cc cross section. The resulting p t differential charm 
cross section is shown in the right panel of Fig. [9] in comparison with the upper and lower limits 
of the corresponding cross section as predicted with a FONLL pQCD calculation 81,82 . The 


upper limit of the calculation is consistent with the measured charm hadron cross sections over 
the full p t range. The same was observed at higher energy with the CDF experiment at the 
Tevatron (see Fig. [6j) and with the ALICE experiment at the LHC (see Fig. 29). 

The pt differential cross section was parametrized with a power-law function (indicated in 
the right panel of Fig. [9]) and extrapolated to zero p t . Furthermore, the rapidity density of 
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the cross section at midrapidity was extrapolated to the full phase space based on calculations 
taking into account the rapidity dependence of charm production, resulting in an extrapolation 
factor of 4.7 ± 0.7. The rapidity density of the cross section and the total cc production cross 
sections were 


^ a cc\ y=0 = 170 ± 45(stat.)^5g(sys.)/ib 


dy 


-59 V 
\ +208 / 


cr C c = 797 ± 210(stat.)_ 29 5(sys.) /ib. 


(3) 

(4) 


3.2.3 Semielectronic heavy-flavor hadron decays 

With both the PHENIX and STAR experiments extensive measurements of electrons from open 
heavy-flavor hadron decays have been conducted at RHIC. The layout of the PHENIX apparatus 
was optimized such that the material budget in the detector aperture was kept at a minimum. 
The measured ratio 155 of electrons of “non-photonic” origin, i.e. mainly electrons from 
open heavy-flavor hadron decays, to electrons of “photonic” origin, i.e. electrons from photon 
conversions in material in the detector aperture and from light meson Dalitz decays, is shown 
in the left panel of Fig. 10 for pp collisions at \fs = 200 GeV. Two complementary techniques 


to separate non-photonic from photonic electrons, i.e. the cocktail subtraction method and the 
converter subtraction method as described in more detail above, gave consistent results for the 
non-photonic/photonic electron yield ratio. For p t above ~ 2 GeV/c this ratio exceeded one. 
The original layout of the STAR apparatus included the rather thick (in terms of radiation 
lengths) SVT detector (see above) which was removed from the experimental setup prior to the 


2008 run of RHIC. As demonstrated in the right panel of Fig. 10 the non-photonic/photonic 
electron yield ratio measured in pp collisions at \fs = 200 GeV with the STAR experiment 
was similar in magnitude to the ratio measured with PHENIX when the SVT was not part 
of the setup 174 . With the SVT in place, this ratio was substantially smaller, rendering the 


measurement of electrons from open heavy-flavor hadron decays very challenging, in particular 
at low p t . 

After the first measurements of electrons from open heavy-flavor hadron decays in pp col¬ 
lisions with the PHENIX and STAR experiments at RHIC, over the following few years the 
precision of the data was improved significantly and the transverse momentum range covered 
was extended both to lower and higher p t . The current status of these measurements is doc¬ 


umented in Fig. [TTJ which shows the invariant p t differential cross sections of electrons from 
heavy-flavor decays measured with the PHENIX experiment (left pan el) [155 in the range 
0.3 < p t < 9.0 GeV/c and with the STAR experiment (right panel) |174|. For the STAR 
experiment two data sets are shown. In the 2005 RHIC run, with the SVT in place, electrons 
were measured in the range 1.2 < p t < 10 GeV/c (triangle symbols). With additional data 
from the 2008 run, in which the SVT was removed from the STAR setup, the precision of the 
measurement could be significantly improved in the range 3 < p t < 10 GeV/c (circle sym¬ 
bols). Corresponding predictions of the electron cross section from heavy-flavor decays from a 


FONLL pQCD calculation 81,82 are compared with the data in the upper panels of Fig. 11 


The ratios of the data and the FONLL calculation are shown in the lower panels of Fig. [TT]on a 
linear scale to allow a better quantitative comparison. The PHENIX and STAR measurements 
are consistent with each other and, within the experimental and the theoretical uncertainties, 
which are dominated by the uncertainties of the chosen normalization and factorization scales, 
the pQCD calculation is in agreement with the data. However, a clear tendency is observed 
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Figure 10: Ratios of electrons from “non-photonic” sources to electrons from “photonic” sources 
measured in pp collisions at yfs = 200 GeV with the PHENIX experiment (left panel) and the 
STAR experiment (right panel) at RHIC. The “non-photonic” electron spectrum is dominated 
by electrons from open heavy-flavor hadron decays. Contributions from other sources do not 
play an important role as demonstrated for the PHENIX case in the left panel (taken from 
Ref. 11551). The important impact of the material budget is shown in the right panel for the 
STAR experiment, in which the material budget was reduced significantly in the 2008 RHIC 
run compared to the 2005 run (see text for details) (reprinted with permission from Ref. 1174 ; 
Copyright (2011) by the American Physical Society). 


that the measured cross sections are close to the upper limit of the predicted cross section, 
consistent with the STAR D-meson measurement (see Fig. [9j). 


According to the FONLL calculation the electron spectrum from heavy-flavor hadron decays 
is completely dominated by charm decays at low p t , which is expected because of the much 
smaller total bottom production cross section compared to charm and, in addition, because of 
the harder bottom decay electron spectrum compared to electrons from charm decays. There¬ 
fore, from the measured p t differential cross section of electrons from heavy-flavor hadron decays 
the total charm production cross section could be determined. The necessary extrapolation of 
the measured electron cross section from the p t and rapidity intervals covered with the PHENIX 
experiment to the full phase space was done in three steps. First, the measured invariant cross 
section was extrapolated to zero p t to obtain the integrated electron cross section per unit 
rapidity at midrapidity over all p t . The p t range covered by the PHENIX data included more 
than 50% of the total electron yield from charm decays at midrapidity. Second, the integrated 
electron cross section was converted into a charm cross section, taking into account the effective 
charm to electron branching ratio and a kinematic correction factor to account for the differ¬ 
ent rapidity distributions of D mesons and their decay electrons. Third, the resulting rapidity 
density of the cross section at midrapidity was extrapolated from midrapidity to the full phase 
space. A systematic uncertainty of 15% was quoted for the extrapolation factor which was 

The resulting rapidity 


6.6. Details concerning this procedure can be found elsewhere 155 


density of the charm production cross section at midrapidity and the total charm production 
cross section from the electron measurement with PHENIX are 
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Figure 11: Invariant p t differential production cross section of electrons from open heavy-flavor 
hadron decays measured in pp collisions at = 200 GeV with the PHENIX experiment (left 
panels) and the STAR experiment (right panels) at RHIC. The upper panels in both plots 
include a comparison of the data with corresponding cross sections from a FONLL pQCD cal¬ 
culation. The lower panels show the ratio of the data and the calculation, where the upper 
(lower) curves indicate the upper (lower) limit of the FONLL calculation. More details are pro¬ 
vided in the text (reprinted with permission from Refs. 1155] (left) and 174| (right); Copyright 
(2011) by the American Physical Society). 


119 ± 12(stat.) ± 38(sys.) pb (5) 

551 ± 57(stat.) ± 195(sys.) pb. (6) 

Unfortunately, the electron data measured with the STAR experiment did not extend to 
low enough electron p t , such that an extrapolation to the full phase space suffered from large 
systematic uncertainties and a sensible total charm production cross section from the electron 
measurement alone could not be quoted. 

With increasing electron p t the contribution from bottom decays grows. According to the 
FONLL calculation, at about 4 GeV/c the contributions from charm and bottom decays are 
similar, and towards higher p t bottom decays become the dominant source of electrons from 
heavy-flavor hadron decays. 

It is important to disentangle the charm and bottom decay contributions experimentally 
such that the FONLL pQCD predictions can be tested separately, and to provide individual 
references for the decay electron spectra from the two heavy quark species in Au+Au collisions 
at RHIC. Both in PHENIX and STAR such separation techniques were developed based on 
electron-hadron correlation measurements. Due to the fact that the combinatorial background 
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Figure 12: Statistical separation of the charm and bottom decay contributions to the electron 
sample from heavy-flavor hadron decays in pp collisions at \fs = 200 GeV with the PHENIX 
experiment (upper two panels) and with the STAR experiment (lower two panels). The invari¬ 
ant mass distribution of electron-hadron pairs measured with PHENIX is compared to charm 
and bottom decay templates obtained from simulations using PYTHIA and EVTGEN for two 
electron p t ranges (upper left panel). The azimuthal angle between electrons from heavy-flavor 
decays and hadrons measured with STAR is compared to charm and bottom decay templates 
obtained with PYTHIA simulations for two electron p t ranges (lower left panel). Relative con¬ 
tributions of bottom decays to the total yield of electrons from heavy-flavor hadron decays as 
functions of the electron p t are compared to predictions from FONLL pQCD calculations (up¬ 
per and lower left panel) (reprinted with permission from Refs. |163j (upper) and |175| (lower); 
Copyright (2009, 2010) by the American Physical Society). 
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in two-particle correlation measurements increases dramatically with particle multiplicity, the 
statistical separation of electrons from charm and bottom decays was possible only in pp colli¬ 
sions at RHIC up to now. 

In the PHENIX case, the extraction of the relative contributions from charm and bottom 
decays was based on the partial reconstruction of the decay D —» e + K~X (and charge conju¬ 
gate) 11631. The invariant-mass distribution of unlike charge-sign electron-hadron pairs exhibits 
a correlated signal below the D-meson mass because of the charge correlation in D-meson de¬ 
cays. The upper left panel of Fig. [I2| shows this invariant-mass distribution normalized to the 
total yield of electrons from heavy-flavor decays in two hadron p t bins. In this analysis, kaon 
identification was not performed but the mass of all reconstructed hadrons was assumed to 
be the kaon mass. Using a combination of the PYTHIA 120 and EVTGEN 1176 simulation 
packages template distributions were generated that described the expected electron-hadron 
invariant-mass distributions from charm and bottom decays, shown as long-dashed and short- 
dashed lines, respectively, in the upper left panel of Fig. [12} A linear combination of these 
template distributions was fitted to the data with the relative contributions to the total dis¬ 
tributions as the only fit parameter. The resulting contribution of bottom decays to the total 


electron spectrum from heavy-flavor hadron decays is shown in the upper right panel of Fig. 02 
as a function of the electron p t . Below 3 GeV/c and above 5 GeV/c only upper and lower 
limits, respectively, of the relative bottom decay contribution could be measured. Nevertheless, 
a clear trend is observed that the bottom decay contribution increases with increasing electron 
p t . Moreover, the measured relative contributions are in good agreement with the expectation 


from a FONLL pQCD calculation |81,82 


In the STAR case, the relative contributions from charm and bottom decays were separated 
statistically from each other via the measurement of angular correlations 175 . In a first 


approach, the distribution of the azimuthal angle A0 between electrons from heavy-flavor 
decays and charged hadrons was studied. Such distributions, normalized to the number of 
electrons from heavy-flavor hadron decays, are shown for two electron p t bins in the lower left 
panel of Fig. [12} The shapes of these distributions, in particular the width of the near side peak 
at A0 = 0 is largely given by the decay kinematics of the charm and bottom hadron decays, 
where the larger mass of bottom hadrons leads to a broader peak compared to charm hadron 
decays. Template distributions reflecting the expected shapes of the electron-hadron azimuthal 
distributions as shown in the lower left panel of Fig. 12 were generated with PYTHIA 120|. The 
contribution of bottom decays relative to the sum of charm and bottom decays was obtained 
from a fit of the measured shape of the angular distribution using a linear combination of 
the charm and bottom templates with the relative normalization as fit parameter. A second, 
independent measurement of the contribution from bottom decays was conducted using the 
angular correlation of electrons from heavy-flavor hadron decays and D° mesons reconstructed 
via their K7T decay. Again, template distributions obtained via a PYTHIA [T20| simulation 
were used to quantify the charm and bottom decay contributions. The observed away-side 
correlation signal in the electron-D 0 azimuthal angle distribution was attributed to prompt 
charm hadron pair production (~ 75%) and bottom hadron decays (~ 25%), while the near 
side correlation signal was mainly due to bottom hadron decays only. The contributions of 
bottom decays to the total electron spectrum from heavy-flavor hadron decays as obtained 
with these two methods are shown in the lower right panel of Fig. [12] as functions of the 
electron p t . Because of the much larger solid angle coverage of the STAR apparatus compared 
to PHENIX, pair correlations could be measured with a much larger geometrical acceptance 
which is a clear advantage for such analyses. Results from the two techniques used in the STAR 
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Figure 13: Invariant p t differential production cross sections of electrons from charm and bottom 
hadron decays, respectively, as measured with the PHENIX (left panels) and STAR experiments 
(right panels) in pp collisions at y/s = 200 GeV. Predictions from a FONLL pQCD calculation 
(shown as lines) are compared with the data on an absolute scale (upper panels) and in form of 


(right); Copyright (2009, 2010) by the American Physical Society). 


the ratio data/FONLL (lower panels) (reprinted with permission from Refs. [1631 (left) and 175 


experiment agree with each other and with prediction from a FONLL pQCD calculation 81,82 
as demonstrated in the lower right panel of Fig. [12| 


From the measured invariant cross section of electrons from heavy-flavor hadron decays at 
y/s = 200 GeV (see Fig. 0 and the relative contributions from charm and bottom decays (see 
Fig. 12) it is straight forward to determine the electron production cross sections from charm 
and bottom hadron decays separately and compare them with the corresponding predictions 
from a FONLL pQCD calculation 181,821. The results from the PHENIX (STAR) experiments 
are shown in the left (right) panel of Fig. 13 For both cases FONLL is compared with the 


data on an absolute scale in the top panels while the ratio of the data to the calculation is 
shown in the lower panels. Given that FONLL provided a decent prediction for both the cross 
section of electrons from all heavy-flavor hadron decays and the for the relative contributions 
of charm and bottom decays to these cross sections it is not surprising that a good agreement 
between the FONLL predictions and the individual electron spectra from charm and bottom 
decays, respectively, was observed both with the PHENIX and STAR experiments. 

Similar to the charm case (see above), the electron production cross sections from bottom 
decays were extrapolated to the full phase to determine the bottom production cross section. 
The extrapolation done for the STAR data was based on the PYTHIA generator 120 in the 
minimum bias mode for pp collisions and in an alternative, exclusive bottom production mode 
in which only leading order processes were included (see Ref. |175| for details). The total bottom 
production cross sections from these two approaches were a bb = 1.34 fib and a bb = 1.83 //!>, 
Corresponding extrapolations of the PHENIX data 1163 using PYTHIA 


respectively 175 
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Figure 14: Electron-positron pair mass distribution from correlated decays of heavy-flavor 
hadrons measured in pp collisions at y/s = 200 GeV with the PHENIX experiment. Arrows 
indicate upper limits at 95% confidence level where the systematic uncertainties were larger 
than the measured pair yields after background subtraction from other sources. Estimated 
contributions from charm decays (thick dashed line), bottom decays (thin solid line), pairs 
from the Drell-Yan process (thin dashed line), and the sum of these (thick solid line) are shown 
for comparison. See text for more details (reprinted from Ref. 1164]), 

and FONLL resulted in 


d<7 bb , 
dy ,2/=0 


a bb 


°-92to:3i(stat.)+J;3g(sys.) pb 
3.2lJ;^(stat.)^ g(sys.) ph. 


( 7 ) 

( 8 ) 


Both the measurements with the STAR and PHENIX experiments are in agreement within 
uncertainties with the total bottom production cross section a bb = /.ib predicted with 

FONLL |8lJ[82|. To a large extent, the substantial uncertainties of the measurements were 
driven by the unavailable cross section measurements for electron production from bottom decay 
below an electron p t of 3-4 GeV/c, indicating that future bottom production measurements have 
to include the low p t region in order to provide a sensitive test of pQCD calculations over the 
full phase space. 

In addition, with the PHENIX experiment the production of correlated electron-positron 
pairs from the decay of heavy-flavor hadron pairs from associated production in a hard scatter¬ 
ing process was investigated in pp collisions at y/s = 200 GeV 1164,465". To statistically isolate 
such electron-positron pairs, in a first step the combinatorial background from uncorrelated 
background pairs was subtracted from the measured electron-positron pair mass distribution. 
The resulting signal distribution was corrected for the reconstruction efficiency. A pair cocktail 
from sources other than heavy-flavor hadron decays was calculated, similar to the single-electron 
case, and it was subtracted from the mass distribution of correlated electron-positron pairs. The 
remaining mass spectrum, shown in Fig. [Ml was interpreted in terms of electron-positron pairs 
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from correlated open heavy-flavor hadron decays. To statistically disentangle these contribu¬ 
tions mass distribution templates were calculated with PYTHIA |120| which were then fitted 
to the data with the relative normalization as free parameters after the contribution from the 
Drell-Yan process was subtracted. As sources of systematic uncertainties in addition to the 
experimental uncertainties, the value of PYTHIA parameters such as the intrinsic k t as well as 
the choice of the parton distribution function, and the uncertainty of the relative abundance 
of charm hadron species were considered. Furthermore, the dynamical correlations between 
heavy-flavor quark-antiquark pairs are not well known, which gave rise to additional system¬ 
atic uncertainties. Within this approach, the total charm and bottom production cross sections 
were determined as 


a c5 = 518 ± 47(stat.) ± 135(sys.) ± 190(model) pb (9) 

a b i = 3.9 ± 2.5(stat.)^2(sys.) /rb, (10) 


where the model dependent uncertainty of the extrapolation to the full phase space is 
quoted separately for the charm cross section. Alternatively, the rapidity density of the charm 
production cross section at midrapidity and the total charm production cross section were 
determined after the bottom decay contribution was subtracted based on a bottom cross section 
of 3.7 pb as obtained from a systematic study I 1 77] of the total bottom production cross section 


as a function of \fs in pp collisions with PYTHIA 120 . The corresponding values are 


118.1 ± 8.4(stat.) ± 30.7(sys.) ± 39.5(model) /ib (11) 

544 ± 39(stat.) ± 142(sys.) ± 200(model) /ib. (12) 

In summary, a consistent picture emerged from the measurement of electrons and correlated 
electron-positron pairs from open heavy-flavor hadron decays with the PHENIX and STAR ex¬ 
periments in pp collisions at = 200 GeV at RHIC. The measured p t differential single 
electron invariant cross sections are consistent with corresponding predictions from FONLL 
pQCD calculations. Also, such FONLL predictions are consistent with the total charm and 
bottom production cross sections obtained from single electron and electron-positron pair mea¬ 
surements with PHENIX and STAR at RHIC. These electron data provide a crucial reference 
for corresponding measurements in nucleus-nucleus collisions. 

3.2.4 Semimuonic heavy-flavor hadron decays 

The PHENIX experiment provides the unique opportunity at RHIC to study the rapidity de¬ 
pendence of open heavy-flavor production via the measurement of electrons at midrapidity (see 
above) and muons at forward rapidity from semileptonic decays of heavy-flavor hadrons. For the 
muon measurement it is crucial to subtract background from primary hadrons which penetrate 
all absorber layers and from other muon sources, i.e. from the light hadron decays before they 
reach the first hadron absorber material and from the decay of light hadrons that are produced 
inside the muon tracker volume of PHENIX. The techniques necessary for the background sub¬ 
traction procedure were developed for the first measurement of muon production at forward 
pseudorapidity (1.5 < \r/\ < 1.8) over the transverse momentum range 1 < Pt < 3 GeV/c in pp 
collisions at y'i = 200 GeV [l66]. The background from misidentified punchthrough hadrons 
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p+pMS = 200 GeV 



Figure 15: Charm production cross section as a function of rapidity in pp collisions at y/s = 
200 GeV measured with the PHENIX experiment using semileptonic open heavy-flavor hadron 
decays to electrons at midrapidity (closed square) and to muons at forward rapidity (closed 
circle). A prediction from FONLL pQCD calculation is shown for comparison (reprinted with 
permission from Ref. 1-50 j ; Copyright (2012) by the American Physical Society). 


was significantly larger for positively charged compared to negatively charged muon candidates, 
mainly because of the relatively large nuclear absorption length for positively charged kaons 
in the absorber material. Therefore, only negatively charged muon candidates were considered 
for the further analysis. After the statistical background subtraction a muon excess remained 
which was attributed to semimuonic decays of open heavy-flavor hadrons, mainly D-mesons 
in the covered muon p t range. A comparison of the p t differential invariant cross section of 
the excess muons with corresponding calculations from PYTHIA and FONLL showed that the 
measured cross sections exceeded significantly both predictions jl66). The experimental sys¬ 


tematic uncertainties, however, were huge and could not be quantified very precisely because 
of the lack of statistics available from the data set recorded during the 2001/2002 RHIC run. 

With new data from the 2005 RHIC run the measurement could be improved signifi¬ 
cantly 1561. The production cross section of negatively charged muons from open heavy-flavor 
hadron decays was measured over an extended p t range (1 < p t < 7 GeV/c) for rapidities 
1.4 < \y\ < 1.9. The shape of the p t differential invariant cross section predicted by a FONLL 
pQCD calculation agreed well with the measurement. Therefore, FONLL was used to extrap¬ 
olate the measured spectrum to zero p t and to calculate from the muon cross section a charm 
cross section at forward rapidity: 


-^I(s/)=1.65 = 139 ± 29(stat.)^g(sys.) pb (13) 

The rapidity dependence of the charm production cross section as predicted from a FONLL 
pQCD calculation is compared to the cross section measured via decay electron at midrapidity 
and decay muons at forward rapidity in Fig. [l5j where the combined muon result from both 
muon spectrometers is shown with a closed circle symbol at negative rapidity and reflected to 
forward rapidity (open circle). Within the substantial experimental and theoretical uncertain¬ 
ties the predicted cross sections are in agreement with the data. Clearly, more precise data 
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Table 2: Rapidity densities da^/dy and total cross sections <r q q for the production of heavy 
quark antiquark pairs qq (here, ’q’ stands for either charm, V, or bottom, ’b’) measured in 
various channels in pp collisions at \f$ = 200 GeV at RHIC. In addition, the p t and y ranges 
covered by the measurements are given. Due to the limited p t ranges, the measured yields 
have to be extrapolated to zero and high p t to determine da^/dy, and the measured fractions 
are listed. Errors quoted are the statistical (first values) and total systematic (second values) 
uncertainties. 


channel 

kinematic coverage 

dcr m /dy 

dqq ref. 


in y and p t 

fraction measured 
in covered y range 



cc D°, D* ± X 

\y\ < 1 

0.6 <p t < 2.0 GeV/c (D°) 
2.0 <p t < 6.0 GeV/c (D*) 

170 ± 45±H lib 
|^67% 

797 ± 210^295 ywb 

170 

cc —y e 1 * 1 T X 

\y\ < 0.35 

119 ± 12 ±38 yh 

551 ± 57 ± 195 pb 

155 


0.3 <Pt< 9.0 GeV/c 

?»55% 



cc —» e + e _ + X 

\y\ < 0.35 

118.1 ±8.4 ±50.0 pb 

544 ± 39 ± 245 pb 

164 


Pt > 0.2 GeV/c 

^80% 



cc —¥ e + e _ + X 

\y\ < 0.35 

not quoted 

518 ± 47 ± 233 ^b 

164 


p e t > 0.2 GeV/c 

?«80% 



cc — p ± X 

1.4 < \y\ < 1.9 

139 ± 29ii yb 

not quoted 

156 


1 < p t < 7 GeV/ c 




bb —» e 1 * 1 + X 

\y\ < 1 

4.0 ±0.5 ±1.1 nb 

1.34 - 1.83 pb 

175 


3.0 < p t < 10.0 GeV/c 

^6.5% 



bb e ± + X 

\y\ < 0.35 

0.92±8;|f±8;ii 

3.2 pb 

163 


3.0 < Pt < 5.0 GeV/c 

^5% 



bb —> e + e _ + X 

\y\ < 0.35 

not quoted 

3.9 ± 2.5^2 pb 

164 


Pt > 0.2 GeV/c 

?a99% 




are needed in order to sensibly discuss the shape of the rapidity distribution of charm in pp 
collisions at RHIC. 

3.2.5 Total heavy-flavor production cross sections at RHIC 

As discussed in this section heavy-flavor production has been studied in detail in pp collisions at 
RHIC in various channels. In addition to the spectral shapes of heavy-flavor observables, which 
are well described by NLO pQCD calculations, the rapidity densities da c5 /dy and da b y,/dy and 
the total heavy-flavor production cross sections a c5 and a hb are the most relevant observables. 
The corresponding values measured in various channels are summarized in Tab. [2} 

3.3 Heavy flavor in d+Au collisions 

Heavy-flavor observables are in agreement with predictions from FONLL pQCD calculations 
in pp collisions at RHIC, providing an experimental and theoretical baseline for correspond- 
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Figure 16: Electrons from heavy-flavor hadron decays in various centrality classes of d+Au 
collisions and in pp collisions measured with the PHENIX experiment at yj snn = 200 GeV. The 
inset shows the ratio of electrons from heavy-flavor hadron decays measured via the converter 
and cocktail subtraction techniques, respectively (left panel). Nuclear modification factor B ( \,\ 
of electrons from heavy-flavor hadron decays in the most central (a) and most peripheral (b) 
d+Au centrality classes (right panel). See text for more details (reprinted with permission from 
Ref. [181 ; Copyright (2012) by the American Physical Society). 


ing measurements in collisions of heavy nuclei in which modifications of such observables are 
expected due to the interaction of heavy quarks with the produced hot and dense medium. 
However, in such collisions also cold nuclear matter effects play a role. Initial state effects, such 
as the modification of the parton distribution functions in nuclei (e.g. shadowing or satura¬ 
tion [178]), broadening of momentum distributions due to the scattering of incoming partons, 
Cronin enhancement 1179,180 of particle production at intermediate and high p t , and energy 
loss of partons in cold nuclear matter are expected to leave their footprint on open heavy-flavor 
observables. In collisions of heavy nuclei it is nearly impossible to cleanly disentangle cold 
matter effects from those induced by the presence of a hot and dense medium. Thus, it is nec¬ 
essary to investigate collisions of protons or deuterons with heavy nuclei. In such collisions the 
production of a hot and dense medium is not anticipated and, therefore, cold nuclear matter 
effects become accessible. 

At RHIC, d+Au collisions were investigated first in an exploratory run in 2003. Data from 
this run firmly established 182-185) that for light hadrons in central Au+Au collisions the 


strong suppression at high p t relative to binary scaled pp collisions and the suppression of 
back-to-back azimuthal angular correlations is due to final-state interactions with the produced 
hot and dense medium and is not induced by cold nuclear matter effects. However, the inte- 
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grated luminosity sampled in this run at RHIC was not sufficient for detailed measurements of 
open heavy-flavor observables. With the STAR experiment the production of D° mesons was 
investigated EH in the transverse momentum range 0.1 < p t < 3 GeV/c, and electrons from 
heavy-flavor hadron decays were measured in the range 1 < p t < 10 GeV/c |l73|. Within large 
uncertainties the d+Au data were consistent with data from pp collisions scaled by the number 
of binary collisions to account for the different collision geometry. The nuclear modification 
factor A<ia of electrons from open heavy-flavor hadron decays was observed to be consistent 
with a moderate Cronin enhancement. However, within substantial uncertainties the electron 
data were also consistent with an R^a of one for all p t . Thus, no significant indication for cold 
nuclear matter effects on open heavy-flavor observables could be claimed based on the data 
from the 2003 run at RHIC. 

The 2008 d+Au run at RHIC constituted a major step forward in terms of the statistics 
available for analysis. For example, with the PHENIX experiment the sampled luminosity of 
80 nb -1 from the 2008 run corresponds to 30 times the sampled luminosity from the 2003 run. 
Using this data set the invariant p t distribution of electrons from heavy-flavor hadron decays was 
measured in d+Au collisions at yj snn = 200 GeV in several centrality classes 181 . The collision 


centrality was determined event-by-event using the information from the PHENIX beam-beam 
counters. Within a Glauber MonteCarlo calculation the mean number of participants, (N part ), 
and the mean number of binary collisions, (Won), was determined for each d+Au centrality 
class. Electrons were reconstructed at midrapidity (|?y| < 0.35) in the transverse momentum 
range 0.85 < p t < 8.5 GeV/c. Background from sources other than open heavy-flavor hadron 
decays was statistically subtracted using the cocktail and converter subtraction techniques as 
discussed in section 13.2.31 

The resulting electron spectra from heavy-flavor hadron decays are shown in the left panel 
of Fig. [16] for minimum bias d+Au collisions, for four d+Au centrality classes, and for pp data 
recorded in the 2008 run at RHIC with the identical experimental setup as it was used for 
the d+Au measurement. The lines represent fits to the previously published pp reference as 


discussed in section section |3.2.3[ scaled by the appropriate number of binary collisions N co n 
for each d+Au centrality classes. The pp data from the 2008 run of RHIC agree very well with 
the pp reference. To facilitate a direct comparison of the d+Au electron spectra with the pp 
reference the nuclear modification factor Raa, defined as 


RdA — 


dN^ A dp t 


(Aeon) x dN*dpt 


(14) 


was calculated for each d+Au centrality class. R,ia is shown as a function of p t for the most 
peripheral and the most central d+Au centrality class in the right panel of Fig. [l6j While in 
peripheral collisions R r \A is consistent with one within large systematic uncertainties, in central 
d+Au collisions an enhancement of the yield of electrons from heavy-flavor hadron decays is 
observed up to a transverse momentum of 5 GeV/c. Measurements in pp collisions suggest 
that this p t range is dominated by charm hadron decays (see section 3.2.3). Because the total 


heavy-flavor production cross section is expected to scale with the number of binary collisions, 
this enhancement at intermediate p t can be interpreted in terms of a pronounced Cronin-like 
broadening 1179 of the heavy-flavor hadron spectra while simultaneously keeping the integrated 
yield unchanged. A similar Cronin broadening was observed in the light quark sector in d+Au 
collisions at RHIC, with the enhancement growing with hadron mass 1186 


These electron data provide clear evidence that open heavy-flavor hadrons are subject to 
cold nuclear matter effects at RHIC. Consequently, such effects are expected to be present in 


36 




















the initial state of nucleus-nucleus collisions at RHIC. However, in the latter case the Cronin- 
like enhancement will be convoluted with effects from the produced hot and dense medium 
such as partonic energy loss leading to a suppression of particle yields at high p t . While 
further experimental input from fully reconstructed open heavy-flavor hadrons and from a clean 
separation of charm and bottom observables is necessary to disentangle all relevant effects at 
RHIC, it is important to keep this first direct evidence for cold nuclear matter modifications of 
an open heavy-flavor observable in mind when interpreting data from nucleus-nucleus collisions 
at RHIC. 


3.4 Heavy flavor in Au+Au collisions 

3.4.1 Hadronic heavy-flavor hadron decays 

The hot medium produced in heavy-ion collisions at RHIC is expected to leave its footprint 
on heavy-flavor observables because heavy quarks interact with this medium while propagating 
through it. In this context, detailed direct measurements of heavy-flavor hadrons in exclusive 
hadronic decay channels in nucleus-nucleus and the reference pp system are highly desirable. 
While D° and D* + measurements with the STAR experiment have become available for pp 


collisions at yfs = 200 GeV recently (see section 3.2.2), the situation is different for Au+Au 
collisions at the same available energy per nucleon-nucleon pair. D° mesons have been recon¬ 
structed with the STAR experiment at low p t (p t < 2 GeV/c) in minimum bias (0 — 80% central) 
Au+Au collisions via their Kir decay [171] . Due to the unfavorably small ratio of the D° signal 
to combinatorial background the extracted yield carried a systematic uncertainty of 40 — 50%. 
The direct measurement of invariant p t differential heavy-flavor hadron yields over an extended 
p t range as a function of the nucleus-nucleus collision centrality will only be possible once data 
become available from the upgraded STAR and PHENIX experiments, which include silicon 
vertex spectrometers capable of separating secondary heavy-flavor hadron decay vertices from 
the primary collision vertex. 


3.4.2 Semielectronic heavy-flavor hadron decays 


The first measurement of open heavy-flavor production in nucleus-nucleus collisions 157 


was 


conducted with the PHENIX experiment using data recorded during the very first run at 
RHIC. In Au+Au collisions at -y/SNN = 130 GeV, electrons from semileptonic heavy-flavor 
hadron decays were statistically separated from electron background from other sources (mainly 
Dalitz decays of light mesons and photon conversions in material) using a cocktail subtraction 
technique as described above. The resulting electron p t spectra from heavy-flavor hadron decays 


are shown in the left panel of Fig. 17 for minimum bias (0 — 92% central) and central (0 — 10% 
central) Au+Au collisions in the range 0.5 < p t < 3 GeV/c. A PYTHIA calculation (120 
the electron spectrum from heavy-flavor hadron decays in pp collisions at y/s = 130 GeV was 


of 


scaled with the nuclear overlap integrals Taa for minimum bias and central Au+Au collisions, 
respectively, and is compared to the data in the left panel of Fig. [17J The PYTHIA parameters 
were tuned such that charm production data from fixed target experiments at the SPS and 
at FNAL as well as single electron data from the ISR were reproduced well. Within large 
uncertainties the scaled PYTHIA calculations agree with the measured electron yields, which 
in this rather low p t range predominantly are due to charm hadron decays. This observation 
suggests that the total charm yield scales geometrically with the nuclear overlap integral or, 
equivalently, with the number of binary collisions when going from pp to Au+Au collisions, as 
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Figure 17: p t differential yields of electrons from open heavy-flavor hadron decays measured with 
the PHENIX experiment in minimum bias and central Au+Au collisions at -^/snn = 130 GeV. 
For comparison corresponding PYTHIA calculations are shown (left panel). Electron produc¬ 
tion cross sections from heavy-flavor hadron decays from PHENIX and ISR measurements for 
the Pi ranges 1.0 < p t < 1.4 GeV/c (open symbols) and p t >1.4 GeV/c (closed symbols), and 
the derived total charm production cross section per binary collision from PHENIX and other 
experiments as a function of \fs. The solid lines and shaded band represent corresponding 
calculations from PYTHIA and an NLO pQCD calculation respectively (right panel). See text 
for more details (reprinted with permission from Ref. 157 ; Copyright (2002) by the American 
Physical Society). 


it is expected for hard probes such as heavy quarks. The measured yield of electrons with p t > 
0.8 GeV/c at midrapidity corresponds to about 15% of the yield integrated over the full p t range. 
Using the spectral shapes in p t and y as predicted from PYTHIA the total charm production 


cross section per binary collision was estimated. The latter is shown in the right panel of Fig. [17 
together with corresponding results from ISR measurements 117,187 as a function of i Js. In 
addition, the y^s dependence of rapidity densities per binary collision of heavy-flavor electron 
yields at midrapidity in two electron p t intervals are depicted. For comparison, corresponding 
results from PYTHIA and from a HVQMNR NLO pQCD calculation 


188 of the total charm 


production cross section are shown. Within substantial uncertainties the PYTHIA and NLO 
pQCD calculations are in agreement with the data. 


While these electron data provided a first glimpse of open heavy-flavor physics in nucleus- 
nucleus collisions, the limited precision of the measurement excluded any detailed conclusion 
concerning the centrality dependence of open heavy-flavor yields, not to speak of the p t distri¬ 
butions of electrons from heavy-flavor hadron decays. 


During the 2001 run period of RHIC a higher statistics electron data sample was recorded 
with reduced systematic uncertainties at midrapidity from Au+Au collisions at Q.S'nn = 200 GeV 
with the PHENIX experiment. Electrons from heavy-flavor hadron decays were measured in 
the p t range 0.4 < p t < 5.0 GeV/c, which is sensitive to charm production only. In a first 
step 1158 , these higher quality data allowed the measurement of the centrality dependence of 


the electron yields from heavy-flavor hadron decays in the p t interval 0.4 < p t < 4.0 GeV/c. 
The background from other electron sources was subtracted using the converter subtraction 
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Figure 18: Electron yield (left axis) and cross section (right axis) from heavy-flavor hadron 
decays divided by the number of binary collisions as measured with the PHENIX experiment in 
Au+Au collisions at 200 GeV in the range 0.8 < p t < 4.0 GeV/c as a function of the number of 
binary collisions, i.e. as function of collision centrality. For comparison the corresponding cross 
section from pp collisions is shown in addition (left panel). Nuclear modification factor Raa of 
electrons from heavy-flavor hadron decays in four Au+Au centrality classes at -^/snn = 200 GeV 
as a function of p t . Error bars (brackets) depict the statistical (systematic) errors from the 
Au+Au measurement. Error boxes show the systematic errors from the pp reference, and the 
band around one indicates the relative systematic uncertainties of the nuclear overlap integrals 
7a A- For the most central collisions the neutral pion Raa is shown for comparison (right panel) 
(reprinted with permission from Refs. 158| (left) and 11601 (right); Copyright (2005, 2006) by 
the American Physical Society). 


technique (see section 3.1.1 for details), which is superior to the cocktail subtraction technique 
at low p t where the signal to background ratio is small. The resulting electron yields per binary 
collision are shown in the left panel of Fig. [18] as a function of the number of binary collisions, 
i.e. as a function of the Au+Au collision centrality. Within the experimental uncertainties the 
electron yields, which are a measure of the number of produced cc pairs, scale with the number 
of binary collisions and in good agreement with the corresponding cross section measured in 
pp collisions at the same y/s. 

The total charm yield scales with the number of binary collisions consistent with point- 
like production in a hard process. Final state interactions of charm quarks with the hot and 
dense medium produced in Au+Au collisions at RHIC are not expected to have an effect on 
the total open charm yield. However, such effects will i nfl uence the momentum distribution 
of charm hadrons and, consequently, their decay products. Electron spectra from heavy-flavor 
hadron decays should be sensitive to such medium modifications for p t > 2.5 GeV/c. While 
the converter subtraction technique is well suited for the measurement of the total yield of 
electrons from heavy-flavor hadron decays (i.e. for low p t electrons), it is of advantage to use 
the complementary cocktail subtraction technique for a precision measurement of the spectral 
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shape towards larger electron p t . 

Models utilizing radiative energy loss via induced gluon radiation as the relevant energy 
loss mechanism for heavy quarks in the hot and dense medium produced in Au+Au collisions 
at RHIC predict only a moderate energy loss of charm and, in particular, bottom quarks at 
transverse momenta up to a few GeV/c (189,190|. In any case, due to the dead cone effect 98 
the radiative energy loss of heavy quarks should be smaller than for light quarks. In view 
of this expectation, the discovery 160 of a substantial suppression of the electron yield from 
heavy-flavor hadron decays at intermediate and high p t in central Au+Au collisions relative to 
the binary collision scaled yields measured in pp collisions came as a complete surprise. 

The right panel of Fig. [18] shows the nuclear modification factor Raa as a function of p t 
for electrons from heavy-flavor hadron decays in four Au+Au centrality classes at +^nn = 
200 GeV 160 . The pp reference was taken from the measurement described in Ref. 159 


Due to small signal to background ratio at low p t the electrons yields at low p t suffer from 
large systematic uncertainties which are reflected by the substantial uncertainties of Raa for 
Pt < 1.5 GeV/c. However, for all centrality classes Raa is consistent with one in the low p t 
region, which includes the dominant fraction of the total yields, in agreement with the observed 
scaling of the total charm yield with the number of binary collisions. Towards higher p t a 
suppression of electrons from heavy-flavor hadron decays develops which is most pronounced 
in the most central collisions at the highest p t . For comparison, the nuclear modification factor 
of neutral pions from 1191 is shown for the case of the 10% most central Au+Au collisions. 
For p t > 3 GeV/c the nuclear modification factor of electrons from heavy-flavor hadron decays 
agrees quantitatively with the Raa of neutral pions in central Au+Au collisions at RHIC. This 
discovery presents a major challenge for models attempting to describe the partonic energy loss 
of heavy quarks traversing the medium created at RHIC. 

Detailed model comparisons became possible with higher statistics electron data from 
heavy-flavor hadron decays recorded in Au+Au collisions at %.S'nn = 200 GeV with the 
PHENIX 155,1621 and STAR experiments 1173 in the 2004 run at RHIC. With the STAR 
experiment electrons were measured in this run in the transverse momentum range 1.2 < p t < 
10 GeV/c. Electron background from photonic sources, i.e. light meson Dalitz decays and 
electrons from photon conversions in material, was measured via an invariant mass technique. 
Remaining background from other sources was subtracted using a cocktail approach. The left 
panel of Fig. [19] depicts the resulting heavy-flavor electron spectra measured with the STAR ex¬ 
periment in three Au+Au centrality classes, d+Au collisions, and the reference pp system 1173 


For comparison results from FONLL pQCD calculations |81|[82 scaled to the pp measurement 
are shown, which indicate that up to 3 GeV/c electrons originate mainly from charm hadron 
decays while towards higher p t bottom hadron decays become increasingly more important. 
With the PHENIX experiment electrons have been measured in the transverse momentum 
range 0.3 < p t < 9.0 GeV/c at midrapidity (|y < 0.35). Electron background from other 
sources was subtracted statistically using the converter method for p t < 1.6 GeV/c and the 
cocktail method for higher p t . The resulting electron spectra from heavy-flavor hadron decays 
for various Au+Au centrality classes 1155,162 are shown in the right panel of Fig. 19 together 
with the corresponding measurement from pp collisions [l61|. The solid lines depict the line 


shape from a FONLL calculation 81,82 normalized to the pp data and scaled with the number 
of binary collisions for each Au+Au centrality class. 

While the binary collision scaled pp reference agrees with the electron yields from heavy- 
flavor hadron decays at low p t for all Au+Au centrality classes this is different at high p t where 
the heavy-flavor electron yields in Au+Au are suppressed relative to the scaled pp reference. 
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Figure 19: Electron p t spectra from heavy-flavor hadron decays measured with the STAR 
experiment in pp, d+Au, and various Au+Au centrality classes at ^/snn = 200 GeV. FONLL 
predictions 181182 are compared with the data for the pp case, for which the right axis indicates 
the corresponding cross section values (left panel). Invariant p t differential distributions of 
electrons from heavy-flavor hadron decays measured with the PHENIX experiment for various 
Au+Au centrality classes and for pp collisions. The solid lines are a from a FONLL calculation 
and the inset shows the signal to background ratio for minimum bias Au+Au collisions (right 
panel) (reprinted with permission from Refs. 1173] (left) and |155] (right); Copyright (2007, 
2011 ) by the American Physical Society). 


The suppression is more pronounced in central collisions compared to more peripheral central¬ 
ity classes. This behavior is quantified in the left panel of Fig. [20] which shows the nuclear 
modification factor Raa of electrons from heavy-flavor hadron decays for two p t ranges as a 
function of the number of participating nucleons, N part , in Au+Au collisions 11621. For the 
range p t > 0.3 GeV/c, which includes more than 50% of the total electron yield from heavy- 
flavor hadron decays, Raa is consistent with one, confirming the binary collision scaling of the 
total heavy-flavor hadron yield in Au+Au collisions at RHIC. For the range p t > 3 GeV/c, Raa 
decreases with increasing centrality reaching a value of ~ 0.4 for the 10% most central Au+Au 
collisions at y%NN = 200 GeV. The suppression of electrons from heavy-flavor hadron decays 
can not be compared directly to light hadron suppression in the same p t range. However, elec¬ 
trons from charm hadron decays with p t > 3 GeV /c originate to a large extent from the decay 
of D mesons with p t > 4 GeV /c 1120]. Therefore, the electron Raa above 3 GeV/c is compared 
with the neutral pion Raa above 4 GeV/c [1911. In this intermediate p t range the trend visible 
in the left panel of Fig. [20] indicates a smaller suppression of heavy-flavor hadrons relative to 
light-flavor mesons in the hot and dense medium produced in Au+Au collisions at RHIC. How¬ 
ever, for a firm conclusion it would be necessary to consider also cold nuclear matter effects, 
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Figure 20: Nuclear modification factor Baa of electrons from heavy-flavor hadron decays 
measured with the PHENIX experiment in Au+Au collisions at ^/snn = 200 GeV as a 
function of the collision centrality (quantified by the number of participants, N pgit ) for the 
electron transverse momentum ranges p t > 0.3 GeV/c and p t > 3.0 GeV/c, respectively. 
The boxes at Baa = 1 indicate the common relative uncertainty from the pp reference for 
p t > 0.3 GeV/c (right box) and p t > 3.0 GeV/c (left box). Baa of neutral pions with 
Pt > 4.0 GeV/c |l9l] is shown for comparison (left panel). Dielectron yield in the interme¬ 
diate mass range 1.2 < m ee < 2.8 GeV/c 2 per binary collision collision as a function of A part 
in comparison with two bands corresponding to different estimates of the contribution from 
charm hadron decays. See text for more details (right panel) (reprinted with permission from 
Refs. 11621 (left) and 1165] (right); Copyright (2007, 2010) by the American Physical Society). 


which might be different for light and heavy-flavor hadrons and which have been measured with 
limited precision only up to now 173, 181| . Furthermore, the interpretation of electron spectra 
from heavy-flavor hadron decays in Au+Au collisions is not straight forward because of the 
decay kinematics involved and, most important, because the relative contributions from charm 
and bottom hadron decays can be different from those measured with substantial uncertainties 
in pp collisions at the same energy 1163,175 . 


A similar general picture emerges from the centrality dependence of dielectron yields per 
binary collision measured with the PHENIX experiment 465 
1.2 < m ee < 2.8 GeV/c 2 as shown in the right panel of Fig. 20 


in the intermediate mass range 
This mass region is dominantly 
populated by dielectrons from correlated charm hadron decays. The measured yields are com¬ 
pared with two extreme scenarios concerning the correlation between the original charm quark 
and antiquark. In the first approach, a PYTHIA calculation for pp collisions 120 is simply 


scaled with the number of binary collisions to the Au+Au case, thus keeping the cc correlation 
intact, i.e. unmodified by the interaction with the hot and dense medium in Au+Au collisions. 
In the second approach, the kinematic cc correlation is assumed to be completely destroyed by 
the interaction with the medium, i.e. the electrons and positrons are sampled randomly from 
the measured electron spectra from heavy-flavor hadron decays 162 . This approach gives rise 


to a considerable softening of the dielectron mass spectra resulting in a reduced dielectron yield 
in the intermediate mass region compared to the PYTHIA scaling as indicated in the right 
panel of Fig. [20} In view of the suppression of the electron yields from heavy-flavor hadron 
decays at high p t the observed binary collision scaling of the intermediate mass dielectron yield 
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Figure 21: Nuclear modification factor Raa of electrons from heavy-flavor hadron decays mea¬ 
sured as a function of p t with the STAR experiment in d+Au and in the 5% most central Au+Au 


collisions at ^/snn = 200 GeV. For comparison the corresponding Raa of charged hadrons 192 
is indicated as well as results from various theoretical model calculations. See text for more 
details (reprinted with permission from Ref. |173 ; Copyright (2007) by the American Physical 
Society). 


for all collisions centralities might point to the cancellation of two competing effects. With in¬ 
creasing centrality the energy loss of charm increases, resulting in a softening of the dielectron 
mass spectra and, consequently, in a reduced dielectron yield in the intermediate mass region. 
This might be compensated by an increasing yield of dielectrons from thermal radiation with 
centrality, which was observed previously with the NA60 experiment 127,128 at the CERN 
SPS (see section 2.1). 

In order to address the suppression of electrons from heavy-flavor hadron decays quanti¬ 
tatively and to relate this with the energy loss of heavy quarks traversing the hot and dense 
medium produced in central Au+Au collisions at RHIC, model calculations of the nuclear mod¬ 
ification factor i?AA(Pt) from various theoretical approaches are confronted with electron data 
from the most central Au+Au collisions. 

Fig. 21 shows such a comparison for the case of the STAR experiment 173 . The nuclear 


modification factor of electrons from heavy-flavor hadron decays in d+Au collisions is consistent 
with a moderate Cronin enhancement albeit with large uncertainties. For the 5% most central 
Au+Au collisions predictions of the heavy-flavor electron Raa are shown from five different 
model calculations. Curve I corresponds to radiative energy loss via a few hard scattering 
processes in the medium calculated within the DGLV formalism 1189 , where an initial gluon 
density of dN g /dy = 1000 was assumed, consistent with hadron suppression in the light quark 
sector. Another radiative energy loss calculation, employing multiple soft collisions within 
the BDMPS formalism 1190 is shown as curve II, where the maximum transport coefficient 
q = 14 GeV 2 /fm was used which is still consistent with light hadron suppression. Both radiative 
energy loss calculations predict less suppression than is observed. Adding collisional, i.e. elastic, 
energy loss to a DGLY-based calculation and introducing path length fluctuations 94 leads to 
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Figure 22: Raa of electrons from heavy-flavor hadron decays measured as a function of p t 
with the PHENIX experiment in the 10% most central Au+Au collisions at 200 GeV [l55| in 
comparison with the corresponding neutral pion Raa from Ref. 19l] (upper panel). Elliptic 


flow magnitude v 2 of electrons from heavy-flavor hadron decays 1155 and of neutral pions 193 
in minimum bias Au+Au collisions (lower panel). 


a larger, but still not sufficient suppression as demonstrated by curve III. In a complementary 
approach (curve IV) heavy-quark energy loss is implemented via elastic scattering involving 
charm and bottom resonance excitations in the medium as well as leading order t-channel 
gluon exchange 101 . In particular towards high p t , where the contribution from bottom 


hadron decays is expected to become important heavy-quark energy loss models consistently 
seem to underpredict the observed suppression. However, given the relatively large uncertainties 
of the STAR electron data it is difficult to draw definite conclusions going beyond the general 
statement that in central Au+Au collisions the yield of electrons from heavy-flavor hadron 
decays is strongly suppressed at high p t , consistent with substantial energy loss of heavy quarks 
in the hot and dense medium produced at RHIC. 

The PHENIX electron data from heavy-flavor hadron decays provide a higher precision and, 
therefore, allow a better discrimination between different theoretical model calculations 155 


The nuclear modification factor Raa of electrons from heavy-flavor hadron decays, shown as a 
function of p t for the 10% most central Au+Au collisions in the upper panel of Fig. [22| is larger 
compared to the neutral pion Raa 191 at low and intermediate p t but towards high p t the 


suppression of the two particle species is of similar size. Further constraints are provided by the 
magnitude v 2 of the elliptic flow of electrons from heavy-flavor hadron decays as compared for 
minimum bias Au+Au collisions to the corresponding neutral pion v 2 1931 in the lower panel of 
Fig. 22. Although for p t above 2 GeV/c the 7r° v 2 is larger than the v 2 of electrons from heavy- 
flavor decays the latter is clearly larger than zero up to a p t of ~ 3 GeV/c. The suppression 
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Figure 23: Heavy-flavor electron Raa in central Au+Au collisions and v 2 in minimum bias 


collisions predicted from a radiative energy loss calculation within the BDMPS framework [190 
(left panel) 


Raa in central Au+Au collisions from a DGLV calculation |94| with radiative 
energy loss only (uppermost band) and with additional collisional energy loss (middle band) 
for electrons from charm and bottom hadron decays. The thin dashed curves correspond to 
charm hadron decays only. The lowest band shows the prediction from a collisional dissociation 
model 1021 (right panel) (reprinted with permission from Ref. [155 ; Copyright (2011) by the 
American Physical Society). 


of electrons from heavy-flavor hadron decays at high p t and the simultaneous participation of 
these hadrons in the collective dynamical expansion of the fireball, as indicated by the non¬ 
zero flow of decay electrons, can be interpreted in terms of strong interactions of heavy flavor 
with the hot and dense medium produced at RHIC. By confronting predictions from various 
theoretical model calculations with the electron data one might hope to unravel the nature of 
the mechanisms responsible for the observed heavy-flavor medium modifications and, at the 
same time, approach a quantitative understanding of the properties of hot QCD matter. 

A first question is whether the measured Raa and v 2 of electrons from heavy-flavor hadron 
decays require heavy-flavor interactions with the medium on a partonic level or whether (pre) hadronic 
interactions are sufficient to explain the observed medium effects. The Raa and v 2 of D mesons 
was calculated in the framework of the hadron-string-dynamics transport model (HSD) [l94| 
which implements heavy-flavor interactions with the medium on a (pre) hadronic level. Only 
at low p t the predicted D-meson Raa is consistent with the electron measurement. Towards 
higher p t the calculated D-meson energy loss is not sufficient to account for the small electron 
Raa■ In addition, the substantial v 2 of electrons from heavy-flavor hadron decays can not 
be reproduced in this model, indicating that strong partonic interactions are relevant for the 
observed medium effect on heavy-flavor observables in Au+Au collisions at RHIC. 

Originally it was expected that the dominant energy loss mechanism of heavy quarks in hot 
QCD matter at RHIC would be induced gluon radiation, leading to a reduced suppression of 
heavy quarks with respect to light quarks due to the dead-cone effect. Predictions using the 
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BDMPS framework 190 , in which radiative partonic energy loss is implemented via multiple 


soft interactions of the partons with the medium, are confronted with the PHENIX electron 
data [155] in the left panel of Fig. [23] The relevant parameter concerning the strength of the 
parton-medium interaction in this model is the transport coefficient q. For the curves shown 
in the left panel of Fig. 23 q = 14 GeV 2 /fm was used which is close to the upper limit still 
consistent with the hadron suppression observed in the light-flavor sector. In particular at 
high p t , the model prediction is difficult to reconcile with the data if bottom hadron decays 
contribute significantly to the electron spectra in this region. Since this BDMPS calculation 
does not include any collective effects at low p t the comparison of the predicted v 2 with the 
data is of limited usefulness. It is interesting to note that the predicted u 2 > 0 provides a lower 
limit of v 2 at high p t which is due to the energy loss in the medium in combination with the 
initial geometry of the collision zone. 

An alternative radiate energy loss approach is implemented in the DGLY framework 194] 189 


in which partons interact with the medium via a few hard collisions. The predicted Raa of 
electrons from heavy-flavor hadron decays, shown as the uppermost band in the right panel of 
Fig. [23] is clearly above the measured Raa at high p t , consistent with the BDMPS prediction. 
It was pointed out |99| that for heavy quarks radiative and collisional energy loss via elastic 
scattering are of similar magnitude over a broad kinematic range at RHIC. Adding collisional 
energy loss to the DGLV calculation leads to a reduction of the electron Raa as demonstrated 
by the middle band in the right panel of Fig. [23] which is still above the measured electron Raa 
at high p t . 


Following a different Ansatz 195-197 , collisional energy loss in pQCD was implemented 


using a running coupling constant a e ff and replacing the Debye mass uid by a hard thermal loop 
calculation. An important parameter in this model is a K factor which is a scaling parameter for 
the collisional energy loss. Raa and v 2 predicted from this model are compared with the data in 
Fig. [24} With a K factor of 1.8, the model prediction is in good agreement with the measured 
Raa of electrons from heavy-flavor hadron decays, not only for the most central collisions but for 
all Au+Au collision centralities. Furthermore, K = 1.8 leaves room for additional contributions 
from radiative or other energy loss mechanisms which are not implemented explicitly. For the 
calculation of the elliptic flow magnitude v 2 the interaction time is important in this model, 
where a later decoupling of heavy quarks from the medium in general can produce larger values 
of v 2 . The curve shown in the right panel of Fig. [24] corresponds to the situation in which the 
hadronization of heavy quarks takes place at the end of the mixed phase. 

Another class of models uses a Langevin-based transport approach to describe the propa¬ 
gation of heavy quarks in the hot medium, i.e. these quarks are placed into a thermal medium 
and the quark-medium interaction is implemented by uncorrelated elastic momentum kicks. 
A relevant parameter in this approach is the heavy-quark diffusion coefficient. The predic¬ 
tions from one of these models 100 are compared with the data for two different values of 
the diffusion coefficient in the left panel of Fig. [25] While the general features of the data are 


reproduced by this model, it can not simultaneously describe the measured Raa and v 2 of single 
electrons from heavy-flavor hadron decays with a single value of the heavy-quark diffusion coef¬ 
ficient. In a similar approach, also based on Langevin transport, the elastic scattering of heavy 
quarks in the medium is mediated by the resonant excitation of charm and bottom states in the 
medium 101,198 . Predictions from this model are confronted with the data for two different 
choices of the resonance widths in the right panel of Fig. [25] Both the measured suppression 


and the azimuthal anisotropy of electrons from heavy-flavor hadron decays are reproduced by 
the model calculation. 
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Figure 24: Nuclear modification factor Raa of electrons from heavy-flavor hadron decays in 
central Au+Au collisions (left panel) and corresponding elliptic flow magnitude v 2 in minimum 
bias collisions (right panel) as predicted from a collisional energy loss model 1195-1971 in com¬ 
parison with the PHENIX data (reprinted with permission from Ref. 1155 ; Copyright (2011) 
by the American Physical Society). 


The diffusion constant, which is the driving parameter of both Langevin-based transport 
models introduced above, can be related to the shear viscosity q of the hot QCD matter 
through which the heavy quarks propagate, where the ratio of q to the entropy density s 
is the quantity that determines the damping rate in the system. In that sense the energy 
loss of heavy quark in the medium and the simultaneous participation of the heavy quarks 
in the medium’s collective expansion can be interpreted as a damping of the initial heavy- 
quark dynamics by the hot and dense medium. Arguments based on the uncertainty principle 
exclude zero viscosity for any thermal system 199 . Furthermore, it was shown |43| that specific 


classes of conformal field theories with gravity duals in the anti-de-Sitter space (AdS/CFT 
correspondence) exhibit a ratio q/s = 1/Att. It was conjectured that this value represents a 
bound for all relativistic thermal field theories [431. As it is discussed in detail in Ref. [155,162 


the Raa and v 2 measurements of electrons from heavy-flavor hadron decays suggest a value of 
q/s (1.33 — 2)/47 t, close to the conjectured lower bound. In turn, AdS/CFT-based energy loss 
models 1104| give a reasonable description of the measured Raa of electrons from heavy-flavor 
hadron decays while the measured v 2 is severely underestimated, q/s can not only be inferred 
from measurements sensitive to the coupling of heavy quarks to the medium as discussed here, 
but also from observables related to the bulk of the medium, involving light flavor quarks only. 
Such estimates based on elliptic flow measurements of light hadrons 1200,2011 and measurements 
of fluctuations 12021 and entropy production |203| are supported by detailed viscous relativistic 
hydrodynamic calculations (204,2051 and they are consistent with the q/s value estimated here 
using electrons from heavy-flavor hadron decays. 

It is an interesting question whether heavy-quark bound states can be formed in hot QCD 
matter. The relevant pQCD dynamics of open heavy-flavor production was worked out in 
Ref. [102 , and the GLV framework was extended to include composite qq systems such that 
the medium induced dissociation probability of D and B mesons propagating through the hot 
and dense medium could be calculated. In this scenario the sequential fragmentation of heavy 
quarks and dissociation of heavy-flavor mesons results in an effective energy loss. Quantitatively, 
the formation time and the detailed expansion dynamics of the produced hot QCD matter as 
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Figure 25: Nuclear modification factor Raa of electrons from heavy-flavor hadron decays in cen¬ 
tral Au+Au collisions and corresponding elliptic flow magnitude v 2 in minimum bias collisions 
as predicted from Langevin-based transport models without 11001 (left panel) and with 
(right panel) resonance excitation of heavy-flavor like hadrons in the medium. See text for de¬ 
tails (reprinted with permission from Ref. EH; Copyright (2011) by the American Physical 
Society). 



well as the in-medium formation time of the heavy-flavor hadrons play a role. The resulting 
nuclear modification factor of electrons from heavy-flavor hadron decays, shown as the lowest 
band in the right panel of Fig. [23| is in reasonable agreement with the measurement as are 
result from improved calculations which combine sequential heavy-flavor hadron formation and 
dissociation with charm and bottom energy loss on the parton level (103 . 

In all above model calculations, addressing the production of heavy-flavor hadron decay 
electrons, the same chemical composition of charm and bottom hadrons is assumed for both 
pp and Au+Au collisions. However, in the light-flavor sector an enhancement of baryon yields 
at intermediate p t relative to meson yields has been observed in central Au+Au collisions at 
RHIC 107 . A similar baryon anomaly in the heavy-quark sector could lead to an enhancement 
of A c relative to D mesons. Given the smaller semileptonic decay branching ratio of charm 
baryons compared to mesons and because of a softer p t spectrum of electrons from charm 
baryon decays, Raa < 1 would be a consequence of such a baryon anomaly [108| even without 
any additional heavy-flavor energy loss in the hot QCD matter. Assuming a charm baryon 
to meson ratio close to one a nuclear modification factor of about 0.65 was predicted 109 


for electrons from heavy-flavor hadron decays. Without additional experimental input it will 
be impossible to resolve this issue. Clearly, the reconstruction of displaced secondary decay 
vertices and exclusive measurements of hadronic heavy-flavor decays are highly desirable as 
they would not suffer from such ambiguities. 

The issue of charm versus bottom energy loss suffers from similar ambiguities. The relative 
contribution, re, from bottom hadron decays to the electron spectra from heavy-flavor hadron 
decays has been measured with substantial uncertainties in pp collisions 1163,1751. For Au+Au 
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Figure 26: 90% confidence level contours for Raa of electrons from D-meson and B-meson 
decays for electron p t > 5 GeV/c in comparison with model calculations. See text for more 
details (reprinted with permission from Ref. |175j; Copyright (2010) by the American Physical 
Society). 


collisions, however, experimental constraints are missing. Therefore, the only practical approach 
to separate the nuclear modification factors of electrons from charm hadron decays, i?AA, an d 
from bottom hadron decays, i?AA, i s to combine the measurement of t~b in pp collisions with 
the measurement of Raa °f electrons from all heavy-flavor hadron decays in Au+Au collisions 
via the equation 


Raa = (1 - r B )R' A D A + r B RZ 


AA- 


(15) 


Since above an electron p t of about 5 GeV /c both re, as measured with the STAR experiment 
and -Raa ) as measured in central Au+Au collisions with the PHENIX 
do not show a significant p t dependence, such an attempt to separate 


in pp collisions 
experiment 


175 


155,162 


the contributions from charm and bottom hadrons decays was conducted in this p t region as 
discussed in detail in Ref. [ 175] . The results are summarized in Fig. [26} which shows the most 


probable values for Rj a and R/% obtained from the data as open circles and the corresponding 
90% confidence limits as dashed lines. Clearly the data imply a substantial suppression not 
only of electrons from charm but also from bottom hadron decays. The predictions from three 
models are confronted with the data in Fig. [26} A DGLV calculation with radiative energy 


loss only (model I) |189| predicts a too small suppression for charm and, in particular, bottom 
energy loss. A collisional energy loss model including collisional dissociation, partonic energy 
loss, and cold nuclear matter effects (model II) [l03| is in reasonable agreement with the data 
as is a model in which collisional energy loss is mediated by resonant elastic scattering (model 
III) [101. Not shown are calculations within the AdS/CFT framework |206|,|207j|, which also 
predict small Raa for electrons from charm and bottom hadron decays. Even more important, 
the ratio of the nuclear modification factors of charm and bottom are predicted 208 to be 


significantly different for AdS/CFT calculations and pQCD-based calculations as shown for 
LHC energies in Fig. [27} emphasizing the importance of separate measurements of charm and 
bottom hadron energy loss in the future. 
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Figure 27: Ratio of the nuclear modification factors of charm and bottom calculated within the 
frameworks of AdS/CFT 1208 and WHDG pQCD [941 with a wide range of input parameters 
as a function of p t (reprinted from Ref. 1208)). 


In summary, the measurements of electrons from heavy-flavor hadron decays in pp and 
Au+Au collisions at RHIC as reviewed here constitute the first systematic application of open 
heavy-flavor observables as a tool to probe the properties of the hot QCD matter produced 
in high-energy nucleus-nucleus collisions. It was shown that the total charm yield scales pro¬ 
portional to the number of binary collisions in nucleus-nucleus collisions at RHIC, i.e. charm 
quark-antiquark pairs are almost exclusively produced in point-like hard scattering processes. 
Additional production processes such as thermal production in the hot medium do not play a 
significant role. The produced heavy quarks propagate through the medium where they are 
subject to strong interactions with the medium, which do not change the overall heavy-flavor 
yields but modify the phase-space distributions to a large extent. Electrons from heavy-flavor 
hadron decays are suppressed at high p t , quantitatively similar to what is observed for light 
hadrons. This observation is in stark contrast with the expectation from energy-loss models 
based on radiative energy-loss mechanisms. Collisional energy loss can not be ignored for heavy 
quarks and/or alternative energy-loss mechanism need to be considered. At the same time, the 
measured v 2 of these electrons indicates that charm quarks participate in the collective expan¬ 
sion of the produced hot QCD medium. These observations support the characterization of 
the hot and dense QCD matter formed in Au+Au collisions at RHIC as a strongly interacting 
medium with partonic degrees of freedom which features a minimal ratio of shear viscosity to 
entropy density and, therefore, was dubbed an almost perfect fluid. However, semielectronic 
heavy-flavor hadron decays provide only indirect information. In particular, a precision sep¬ 
aration of charm from bottom is difficult via this approach, and the chemistry of charm and 
bottom hadron states is completely inaccessible. These issue can only be addressed with the 
upgrades of the PHENIX and STAR experiments. 


3.4.3 Semimuonic heavy-flavor hadron decays 

While extensive systematic studies of heavy-flavor production in heavy-ion collisions have been 
conducted at RHIC in the semielectronic decay channel at mid-rapidity, corresponding mea¬ 
surements in the semimuonic decay channel at forward rapidity have become available only 
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Figure 28: p t -differential invariant yields of negatively charged muons from heavy-flavor hadron 
decays in pp and three centrality classes of Cu+Cu collisions measured at forward rapidity 
with the PHENIX experiment at yj snn = 200 GeV. The solid lines depict a fit to the pp data 
scaled by the appropriate number of binary collisions to the Cu+Cu centrality classes (left 
panel). Nuclear modification factor Raa for negative muons from heavy-flavor hadron decays 
measured as a function of p t in the 20% most central Cu+Cu collisions. For comparison, a 
prediction from a collisional dissociation model |103| is shown (middle panel). Comparison of 
the nuclear modification factors of heavy-flavor muons at forward rapidity (1.4 < y < 1.9) and 
p t > 2 GeV/c in Cu+Cu collisions 11561 and heavy-flavor electrons at mid-rapidity (|t 7| < 0.35) 
and p t > 3 GeV/c in Au+Au collisions 1162| as a function of the collision centrality, expressed by 
the number of participating nucleons (right panel) (reprinted with permission from Ref. [156 
Copyright (2012) by the American Physical Society). 


recently. 

With the PHENIX experiment negatively charged muons from the decay of heavy-flavor 
hadrons have been measured in the rapidity range 1.4 < y < 1.9 and for transverse momenta 
1 < pt < 4 GeV/c in three centrality classes of Cu+Cu collisions at %.S'nn = 200 GeV 
Background, mainly from light hadrons, was estimated using a MonteCarlo approach and it 
was subtracted statistically from the inclusive muon candidate sample. The resulting invariant 
yields measured in Cu+Cu collisions are compared as a function of the muon p t with the 
corresponding data from pp collisions in the left panel of Fig.[28| where the solid lines correspond 
to a fit of the pp data scaled with the appropriate numbers of binary collisions to the Cu+Cu 
centrality classes. While the muon yields from heavy-flavor hadron decays in peripheral and 
mid-central Cu+Cu collisions are in reasonable agreement with the binary collision scaled pp 
reference, for the 20% most central Cu+Cu collisions a suppression of such muons is observed 
at high pt relative to the pp fit. 


156 


This suppression is quantified via the nuclear modification factor Raa of muons from heavy- 
flavor hadron decays shown as a function of p t in the middle panel of Fig. [28] for the 20% 
most central Cu+Cu collisions. A prediction from a theoretical model 103 agrees well with 
the measurement at high p t . This specific model includes elastic and inelastic in-medium 
energy loss of heavy quarks, the dissociation of heavy mesons in medium, and cold nuclear 
matter effects such as shadowing and initial state energy loss due to multiple scattering of 
the incoming partons |103|. A related model prediction 11021 is also in reasonable agreement 
with the measured nuclear modification factor of electrons from heavy-flavor hadron decays in 
Au+Au collisions at mid-rapidity (see Fig. 23). 
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As demonstrated in the right panel of Fig. [28j the suppression of muons from heavy-flavor 
hadron decays at forward rapidity in central Cu+Cu collisions is comparable to the suppression 


observed in central Au+Au collisions for electrons from heavy-flavor hadron decays 155,162 


This is a surprising observation since the Bjorken energy density of hot QCD matter produced 
at mid-rapidity in central Au+Au collisions at RHIC is expected to be larger by at least 
a factor of two compared to the energy density in the forward rapidity region in Cu+Cu 
collisions. Consequently, the substantial suppression observed in Cu+Cu collisions might point 
to significant cold nuclear matter effects playing a role in particular at forward rapidity. 

Major progress at forward rapidity requires improved heavy-flavor production measurements 
at RHIC. In particular, the PHENIX forward silicon vertex detector will allow the separation 
of charm and bottom contributions to the muon sample from heavy-flavor hadron decays. 


4 Open heavy-flavor measurements at the LHC 


4.1 Experiments studying heavy-flavor production at the LHC 

The LHC delivered the first pp collisions at injection energy, i.e. at \fs = 0.9 TeV in November 
2009. While data were recorded with all experiments during the commissioning phase at this 
energy, the accumulated statistics was not sufficient to investigate the production of heavy- 
flavor hadrons. Substantial data sets, providing already detailed information on heavy-flavor 
production in pp collisions, were recorded by all four large LHC experiments at an energy of 
y/s = 7 TeV in the years 2010/2011. In particular, these data provided a reference for the first 
Pb+Pb collision data recorded at the LHC at y^NN = 2.76 TeV in November/December 2010 
and 2011 with the ALICE, ATLAS, and CMS experiments. The LHCb experiment does not 
participate in heavy-ion runs. Because of the difference in the available energy per nucleon- 
nucleon pair the 7 TeV pp data have to be scaled to the lower energy of 2.76 TeV to provide 
a reference for measurements in Pb+Pb collisions. As will be demonstrated in Section m 
all heavy-flavor production data from 7 TeV pp collisions at the LHC are in good agreement 
with predictions from FONLL pQCD calculations 83 . Therefore, FONLL driven yfs scaling 


procedures could be used to determine pp references at the proper energy for Pb+Pb collisions. 
In order to validate the FONLL based yfs scaling and to provide experimental references for 
observables which can not be scaled by pQCD model calculations, the LHC delivered pp colli¬ 
sions at the energy of ^Js = 2.76 TeV for a short period of a few days in March 2011. While 
the LHC operates at a collision energy of yfs = 8 TeV since April 2012 relevant results from 
this experimental pp campaign are not yet available. 

This review discusses two issues related to the LHC experiments. First, it is demonstrated 
that perturbative QCD calculations are able to reasonably well describe heavy-flavor hadron 
production in pp collisions at the LHC, using data from all four large experiments in the 
hadronic and semileptonic decay channels at mid-rapidity and at forward rapidity. It is ac¬ 
tually a big success of the experiments as well as of the theoretical calculations that this can 
be concluded already during the very first years of LHC operations. It took years to resolve 
discrepancies that seemingly existed between pQCD calculations and measurements at the Fer- 
milab Tevatron. Also at RHIC, long standing discrepancies between the electron spectra from 
heavy-flavor hadron decays measured with the PHENIX and STAR experiments, respectively, 
could be resolved only after a while. Second, the first published results on open heavy-flavor 
production in Pb+Pb collisions obtained with the ALICE and CMS experiments are summa- 
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rized. While these measurements confirm the substantial suppression of heavy-flavor probes 
at high p t with respect to binary collision scaled pp measurements as observed at RHIC, they 
also represent a major step forward not only in terms of the energy scale. For the first time, 
hot QCD matter effects on exclusive heavy-flavor hadron observables can be investigated. In 
Pb+Pb collisions, D mesons have been reconstructed with the ALICE experiment and with 
the CMS experiment J/^ mesons originating from b-hadron decays have been studied. Such 
measurements have only been possible because all of the four large LHC experiments include 
precision silicon vertex detectors, which allow the experimental separation of displaced sec¬ 
ondary vertices from heavy-flavor hadron decays from the primary collision vertex. Augmented 
by inclusive measurements, e.g. of electrons and muons from semileptonic heavy-flavor hadron 
decays, these first exclusive results are only the tip of the iceberg of near time future results 
to be expected from systematic heavy-flavor production studies with all experiments at the 
LHC. In the following, the experimental approaches of the four large LHC experiments to 
open heavy-flavor measurements are summarized briefly, with special emphasis on the ALICE 
experiment. 


4.1.1 The ALICE approach to open heavy-flavor measurements 


The ALICE experiment 1209 was designed and optimized to cope with the large particle mul¬ 
tiplicities in central Pb+Pb collisions at the LHC. The ALICE apparatus consists of two main 
parts: a central barrel detector system housed inside the L3 solenoid magnet, which provides a 
field of 0.5 T, and a forward muon spectrometer (—4 < rj < —2.5) including a dipole magnet, 
which provides a field integral of 3 Tm. This setup is augmented by several smaller auxiliary 
detectors for event characterization and for triggering purposes. Special emphasis was put on 
excellent tracking and particle identification capabilities, in particular also at low momenta. 
With these capabilities excellent measurements of heavy-flavor production are possible through 
the reconstruction of heavy-flavor hadrons in their hadronic decay channels at mid-rapidity as 
well as via semileptonic decays of charm and bottom hadrons both at mid-rapidity and forward 
rapidity. 

The central barrel tracking system consists of the silicon Inner Tracking System (ITS) 
and the Time Projection Chamber (TPC). In the pseudorapidity range \rj\ < 0.9 tracks are 
reconstructed with a momentum resolution better than 4% for p t < 20 GeV/c. The distance of 
closest approach of tracks to the interaction vertex is measured with a resolution better than 
75/um for p t > 1 GeV/c in the plane transverse to the beam direction 1210 . 


Charged hadrons are identified via their specific energy loss dE/dx in the TPC and a time- 
of-flight (TOF) measurement in the ALICE TOF detector. The same techniques allow a clean 
electron identification in the range p t < 6 GeV/c. At higher p t , the Transition Radiation 
Detector (TRD) and, alternatively, the Electromagnetic Calorimeter (EMCal) are necessary 
for the measurement of a pure electron sample. 

The measurement of D mesons at mid-rapidity (\y\ < 0.5) is based on the selection of 
displaced-vertex topologies 210 . Tracks are selected which originate from a secondary vertex 


consistent with a large decay length as expected from the D-meson lifetimes. In addition, a good 
alignment of the reconstructed D-meson momentum with the flight line from the collision vertex 
to the secondary vertex is required. Since the relevant D-meson decay channels (D° —» K _ 7 t + , 
D+ —> K 7r + 7r + , and D* + —»D° 7r + ) involve a kaon as decay product the identification of the 
charged kaon in the TPC and TOF detectors helps to reduce the background. In an invariant 
mass analysis the raw D-meson yields are determined which, then, are corrected for geometrical 
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acceptance and for reconstruction and particle identification efficiency based on a detailed 
simulation of the apparatus. Feed down from B-meson decays is estimated from a FONLL 
pQCD calculation |83| and it is subtracted. 


Bottom production can be investigated via the measurement of 3/ip mesons from the decay of 
bottom hadrons, c/.Tab.[lJ 3/ip mesons are reconstructed via their e + e _ decays at mid-rapidity 
(\y\ < 0.9). The fraction of the total 3/ip yield originating not from the primary collision vertex 
but from a displaced bottom-hadron decay vertex is determined via the measurement of the 
pseudo-proper decay length of the 3/ip mesons 1211 . 


Electrons from heavy-flavor hadron decays are measured at mid-rapidity 12121 in a two step 
procedure. First, the inclusive electron spectrum is determined. Electron candidate tracks 
are identified with the TPC, TOF, TRD, and EMCal detectors. The small remaining hadron 
background is estimated via fits of the TPC dE/dx distribution in momentum slices and it is 
subtracted from the electron candidate spectra. After corrections for geometrical acceptance 
and tracking and particle identification efficiencies derived from simulations the inclusive elec¬ 
tron spectrum is obtained. From this inclusive spectrum a cocktail of electrons from sources 
other than heavy-flavor decays is subtracted. Similar to corresponding analyses at RHIC, the 
most important components of this background cocktail are Dalitz decays of the tt° and q 
meson as well as photon conversions in the beam pipe and in the first pixel layer of the ITS. 
In addition, decays of light vector mesons (p,u,<p), heavy quarkonia (J/ip, Y), and direct ra¬ 
diation from hard scattering processes are included. The cocktail input is based on the 7r° 
and // spectra measured with the ALICE experiment 1213 , the known material budget of the 
apparatus 1214), m t scaling for the p t spectra and yield of light vector mesons, heavy quarkonia 
measurements with the ALICE 1215 and CMS experiments 216, 217 , and next-to-leading- 
order pQCD calculations 1218, 219| for the electron spectra from direct radiation. With the 
ALICE ITS, electrons originating from displaced secondary bottom-hadron decay vertices can 
be separated from decay vertices of shorter lived charm hadrons and from the primary collision 


vertex. This technique allows the measurement of electrons from bottom-hadron decays 220 


Muons with momenta above 4 GeV/c are identified and their momenta are measured in 
the pseudorapidity range —4 < q < —2.5 with the ALICE forward muon spectrometer, which 
consists of a 10 interaction lengths passive front absorber, a small angle beam shield, a set 
of tracking chambers in front of, in, and behind a 3 Tm dipole magnet, and trigger chambers 
located behind an iron wall with a thickness of 7.2 interaction lengths. From the inclusive muon 
candidate spectrum three main background sources have to be removed: muons from the decay 
of primary light hadrons, muons from the decay of hadrons produced in secondary interactions 
in the front absorber, and hadrons that punch through the front absorber. Contributions from 
the decay of W ± and Z° vector bosons are negligible in the muon p t range 2 < p t <12 GeV/c. 
Punch through hadrons are rejected efficiently requiring that muon candidate tracks match 
a track in the muon trigger system. Decay muons are the dominant source of background 
in the relevant p t region. In pp collisions (221) this contribution can be subtracted using 
information from detailed simulations. In Pb+Pb collisions 222 the presence of unknown 


medium effects prevents such a MonteCarlo approach. Instead, the decay muon contribution is 
estimated by extrapolating the p t spectra of pions and kaons measured in Pb+Pb collisions at 
mid-rapidity 223 to forward rapidity and calculating the corresponding decay muon spectra. 
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4.1.2 The ATLAS approach to open heavy-flavor measurements 

The ATLAS collaboration operates one of the two large general purpose detectors at the LHC 
designed to extend the frontiers of particle physics by exploiting the unprecedented high collision 
energy and luminosity delivered by the LHC. The main goals of ATLAS are related to the search 
for the Standard Model Higgs boson, precision tests of QCD, electroweak interactions, flavor 
physics, and the search for physics beyond the Standard Model in pp collisions. Utilizing mainly 
hard probes such as jets, ATLAS also studies the Pb+Pb collisions provided by the LHC. In 
the following, a brief characterization of those subsystems of the ATLAS apparatus is given 
which are relevant for the results published on open heavy-flavor production in pp collisions up 
to now. ATLAS has not yet published open heavy-flavor data from Pb+Pb collisions. 

The ATLAS apparatus covers almost the full solid angle around the interaction point. An 
inner detector (ID) tracking system consisting of a 3 layer silicon pixel detector, a 4 layer 
double-sided silicon strip detector, and a 73 layer transition radiation straw tube tracker is 
located in a 2 T solenoidal magnetic field, and it covers the pseudorapidity range \rj\ < 2.5. 
The ID is surrounded by a high-granularity liquid-argon electromagnetic calorimeter (EM), 
also covering \rj\ < 2.5, and hadron calorimetry provided by an iron-scintillator tile calorimeter 
in the mid-rapidity region and copper/liquid-argon endcap calorimeters towards forward and 
backward rapidity. The outer muon spectrometer (MS) is located inside a toroidal magnetic 
field providing a bending power of 2.5 Tm in the barrel and 5 Tm in the end caps. Muons 
are tracked in the region |? 7 | < 2.7 using monitored drift tube chambers and higher granularity 
cathode strip chambers in the end caps. Muon trigger chambers employ resistive plate chambers 
for \rj\ < 1.05 and thin gap chambers in the region 1.05 < |?y| < 2.4. This general purpose setup 
provides ample opportunities for the measurement of open heavy-flavor observables. A detailed 
descriptions of the detector layout of the ATLAS experiment can be found elsewhere 1224 . 


4.1.3 The CMS approach to open heavy-flavor measurements 

The CMS collaboration operates the second large general purpose detector at the LHC. The 
main goals of CMS are identical to the ones of ATLAS (see Sec. 4.1.2). Open heavy-flavor 


measurement in pp and Pb+Pb collisions with the CMS experiment at the LHC have used 
the inner tracker and the muon system of the apparatus up to now. Therefore, these detector 
components are briefly characterized below. A detailed descriptions of the detector layout can 
be found elsewhere |225 


The central component of the CMS apparatus is a superconducting solenoid with an in¬ 
ternal diameter of 6 m, providing at 3.8 T magnetic field. Immersed in the magnetic field 
are the inner tracker surrounded by a lead-tungstate crystal electromagnetic calorimeter and 
a brass/scintillator hadron calorimeter. The silicon pixel tracker comprises three barrel layers 
and two endcap disks at each barrel end. This device is surrounded by a silicon strip tracker 
consisting of 10 barrel layers and 12 endcap disks at each barrel end. Muons are detected in the 
pseudorapidity interval \r/\ < 2.4 with three types of gas-ionization detectors embedded in the 
magnet’s steel return yoke: drift tubes, cathode strip chambers, and resistive plate chambers. 
Due to the strong magnetic field and the silicon tracker’s high granularity a p t resolution of 
about 1% is reached for muon candidates which are matched with reconstructed tracks. The 
silicon tracker also provides the position of the primary collision vertex with an accuracy of 
about 20 /rm allowing the separation of primary tracks from those originating from displaced 
secondary decay vertices. Weak decays of hadrons, e.g. Kg —» 7r + 7r - , can be reconstructed 
via an invariant mass analysis of charged particle tracks pointing back to a common decay 
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vertex and fulfilling decay topology cuts, even without the identification of the decay hadron 
species. In addition, the CMS experiment includes substantial forward calorimetry in the pseu¬ 
dorapidity range 2.9 < r/ < 5.2 as well as auxiliary forward detectors for triggering purposes 
and beam-halo rejection. In Pb+Pb collisions information provided by the forward detectors 
is used for event selection and the event-by-event centrality determination. As in the case of 
ATLAS, the CMS experiment gives access to open heavy-flavor observables in a multitude of 
channels, not all of which have been addressed in publications yet. 


4.1.4 The LHCb approach to open heavy-flavor measurements 


The LHCb experiment was designed as a single-arm forward spectrometer 12261 aiming at 


the measurement of CP violating and rare decays of heavy-flavor hadrons in pp collisions at 
the LHC. Pb+Pb collisions are not investigated with the LHCb experiment. Nevertheless, in 
the context of this review LHCb provides an important testing ground for pQCD calculations 
regarding the rapidity dependence of heavy-flavor production cross sections in pp collisions at 
the LHC. 


The LHCb apparatus covers the pseudorapidity range 2 < r] < 5. The detector includes 
a high precision tracking system comprising a silicon strip vertex detector surrounding the 
nominal pp interaction point, a large-area silicon strip detector located upstream of a dipole 
magnet providing a bending power of 4 Tm, and three stations of silicon strip detectors and 
straw drift tubes located downstream. With this tracking system a momentum resolution A p/p 
is achieved between 0.4% at 5 GeV/c and 0.6% at 100 GeV/c. For high p t tracks, the impact 
parameter is measured with a resolution of 20 pm. Two ring imaging Cherenkov detectors 
are used for hadron identification. Photon, electron, and hadron candidates are selected with 
a calorimeter system including scintillating pad and pre shower detectors, an electromagnetic 
and a hadronic calorimeter. Alternating layers of iron absorbers and multiwire proportional 
chambers compose the muon detection system. 

In addition to a variety of results concerning CP violating and rare heavy-flavor hadron de¬ 
cays, the LHCb collaboration has also published baseline cross sections of bottom production in 
pp collisions at Qi = 7 TeV, which can be compared with predictions from pQCD calculations. 


4.2 Reference measurements in pp collisions 

As for the RHIC case, open heavy-flavor measurements in pp collisions at \fs = 7 TeV at the 
LHC provide a crucial testing ground for higher order pQCD calculations in a yet unexplored 
high energy regime and, at the same time, they deliver a baseline for corresponding measure¬ 
ments in Pb+Pb collisions, where the interaction of heavy-flavor particles with the produced 
hot and dense QCD medium is expected to leave its footprint on open heavy-flavor observables. 


4.2.1 Hadronic heavy-flavor hadron decays 


The only measurement of mid-rapidity charm hadron production cross sections in pp collisions 
at the LHC via the full reconstruction of hadronic decay channels has been conducted with the 

D°, D + , and D* + /^-differential production cross sections measured in 

Theoretical calculations based on higher-order 


ALICE experiment |210 


the range \y\ < 0.5 (210 


are shown in Fig. 29 


pQCD calculations (FONLL |83| and GM-VFNS 84]) are in reasonable agreement with the 
data within substantial experimental and theoretical systematic uncertainties. 
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Figure 29: p t -differential production cross sections of D° (left panel), D + (middle panel), and 
D* + (right panel) measured at mid-rapidity (|y < 0.5 ) with the ALICE experiment in pp col¬ 
lisions at t/s = 7 TeV |210| in comparison with pQCD calculations [83,84]. The lower panels of 
the individual plots show the ratio of the measured cross sections and the corresponding predic¬ 


tions from higher order FONLL and GM-VFNS pQCD calculations (reprinted from Ref. 210 
with kind permission from Springer Science and Business Media). 





Figure 30: p t - (upper panels) and rapidity-differential (lower panels) production cross sections 
of B° (left panels), B + (middle panels), and B° (right panels) mesons measured with the CMS 
experiment in pp collisions at s = 7 TeV. Predictions from higher order pQCD calculations 
(MC@NLO |230|) are compared with the data (reprinted with permission from Refs. 1227-229 
Copyright (2011) by the American Physical Society). 
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The situation is similar for the production of bottom hadrons in pp collisions at the LHC as 


demonstrated in Fig. 30 which shows the production cross sections of B°, B + , and B° mesons as 
a function of p t and, for high p t mesons, as a function of rapidity [227-2291. Calculations within 


the framework of the higher-order MC@NLO pQCD model [230 are in reasonable agreement 
with the data although both the experimental and the theoretical systematic uncertainties are 
substantial again. 


B + meson production in pp collisions was also investigated with the LHCb experiment 231 


in the kinematic range 0 < p t < 40 GeV/c and 2.0 < y < 4.5. The measured B + cross sections 
are in good agreement with FONLL pQCD calculations 83 . Furthermore, the pseudorapidity 


dependence of bottom hadron production cross sections was measured with LHCb 12321 in the 
range 2 < y < 6 without an exclusive identification of the various hadron species. Also in this 
case, higher order pQCD calculations in the FONLL |83| and MCFM 12331 frameworks are able 
to describe the measured data. 

An alternative approach to bottom hadron production is the measurement of non-prompt 
J /if) mesons originating from displaced bottom hadron decay vertices. All four large experiments 
at the LHC include vertex spectrometers with sufficient resolution for such measurements. The 
relative contribution / B of bottom hadron decays to the inclusive J/A yield increases at mid¬ 
rapidity from ?»10% at low J/c pt to ~2/3 at high p t as demonstrated in the upper left panel 
of Fig. 31 showing data from the ALICE |21l], ATLAS [234 , and CMS [235 experiments 
measured in pp collisions at yfs = 7 TeV. The dashed line superimposed to the data indicates 
a semi-phenomenological function which is used to extrapolate the measured fraction down 
to zero p t - As an example, the ATLAS measurement of the p t -differential production cross 
section on non-prompt J/0 mesons |234| is shown for four different rapidity intervals between 
mid-rapidity and 2.0 < y^/^ < 2.4 in the lower four panels of Fig. 31, Only in the rapidity 
interval 1.5 < yj/$ < 2.0 the ATLAS measurement reaches down to a J/A meson p t of 1 GeV/c, 
whereas the minimum p t in the other rapidity intervals is 5 GeV/c or even higher. FONLL 
predictions are in agreement with the measured /^-differential production cross section for all 
rapidity selections. Furthermore, the measured rapidity dependence of the J/^ cross section 
from bottom hadron decays is reproduced by FONLL within experimental and theoretical 


uncertainties as is demonstrated in the upper right panel of Fig. 31, which includes data from 
ALICE [211 , CMS [235], and LHCb 12361. ALICE has measured the p t -differential cross section 
of non-prompt J /ip meson production in the transverse momentum range p t > 1.3 GeV/c 
at mid-rapidity ([?/[ < 0.9). The shape of the corresponding FONLL prediction was used 
to extrapolate the measurement to zero p t - A corresponding procedure was applied for the 
CMS measurement in the range 1.2 < \y\ < 1.6. The systematic uncertainties related to this 
extrapolation procedure are indicated by the slashed boxes in the upper right panel of Fig. [31] 
Neither the CMS data in the interval 1.6 < \y\ < 2.4 nor the LHCb data had to be extrapolated 
since non-prompt J/A mesons have been measured down to zero p t in these cases. The cross 
section measured with the LHCb experiment at forward rapidity (closed symbols) have been 
reflected around y = 0 (open symbols). The corresponding FONLL calculation |83| and its 
uncertainty are indicated by the solid and dashed lines, respectively. 


4.2.2 Semielectronic heavy-flavor hadron decays 

Both ALICE and ATLAS have published first cross section measurements of electrons from 
semileptonic decays of heavy-flavor hadrons in pp collisions at ^s = 7 TeV at the LHC. The 


ALICE measurement 212 is restricted to a narrow interval around mid-rapidity (|y| < 0.5). 
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Figure 31: Fraction of J/^ from the decay of bottom hadrons in the inclusive J/'ip sample as 
measured with ALICE [211] , ATLAS (234 , and CMS [235 in 7 TeV pp collisions at mid-rapidity. 
Rapidity-differential production cross section of J /ip from bottom-hadron decays measured with 
ALICE 12111, CMS 1235 , and LHCb 12361 (upper two panels, reprinted from Ref. 12111 with 


kind permission from Springer Science and Business Media). p t -differential non-prompt J/'ip 
production cross sections in four rapidity intervals as measured with the ATLAS experiment 
(lower four panels, reprinted from Ref. 1234)). FONLL pQCD predictions of J/L cross sections 
from bottom-hadron decays are compared with the measured differential J/xp cross sections in 
all but the upper left panel, (see text for further details.) 
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Figure 32: p t -differential production cross sections of electrons from heavy-flavor hadron decays 
measured around mid-rapidity with the ALICE 12121 and ATLAS experiments 1237 
collisions at \fs = 7 TeV (left panel, reprinted with permission from Ref. 


212 


m pp 

Copyright 

(2012) by the American Physical Society). /^-differential production cross sections of electrons 
from charm- and from bottom-hadron decays measured with ALICE at mid-rapidity in pp 
collisions at \/s = 7 TeV (right panel, reprinted from Ref. 12201). In both panels corresponding 
predictions from a FONLL pQCD calculation 83 are shown for comparison. 


With ATLAS, electrons were reconstructed |237| in the larger range |y < 2, excluding the 
interval 1.37 < \y\ < 1.52 in which the electron identification and the resolution of the energy 
measurement is degraded by a large amount of material in front of the first detector elements. 
The resulting/^-differential production cross sections, normalized per unit rapidity, are depicted 
in the left panel of Fig. [32} While the ALICE data include most of the total electron cross 
section from heavy-flavor hadron decays, the ATLAS data extend the measurement to higher 
p t . Exploiting the excellent vertexing capabilities of the ALICE ITS, electrons from charm and 
bottom hadron decays can be separated from each other in pp collisions at \fs = 7 TeV 12201. 
The resulting /^-differential production cross sections are shown in the right panel of Fig. |32| 
Predictions from a FONLL pQCD calculation [83| are in good agreement both with the inclusive 
electrons cross sections from heavy-flavor hadron decays as well as with the individual cross 
sections of electrons from charm- and bottom-hadron decays, respectively, within experimental 
and theoretical uncertainties. 
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Figure 33: p t -differential production cross sections of leptons from heavy-flavor hadron decays 
measured with the ATLAS experiment in pp collisions at \/S = 7 TeV at the LHC. Electron 
and muon production cross sections for \p\ < 2.0, excluding the region 1.37 < r/1 < 1.52 (left). 
Muon production cross section for \r/\ < 2.5 (right). Predictions from FONLL and other pQCD 
calculations are compared with the data. Ratios of the measured cross sections and of other 
predictions to the FONLL calculation are shown in the lower panels (see text for further details) 
(reprinted from Ref. HD). 


4.2.3 Semimuonic heavy-flavor hadron decays 


Both with the ALICE |22l] and ATLAS |237| experiments not only electrons but also muons 
from open heavy-flavor hadron decays have been measured at the LHC. The left panel of 


Fig. 33 shows a comparison of the muon p t -differential production cross section from ATLAS 


with the electron data (also shown in Fig. 32) for p t > 7 GeV/c in the pseudorapidity range 
\i]\ < 2.0, excluding the interval 1.37 < | r]\ < 1.52 in which no electron measurement is 
available. In this high p t range the mass difference between electrons and muons is irrelevant. 
Therefore, one would expect the cross sections of electrons and muons from semileptonic heavy- 
flavor hadron decays to be equal which is confirmed within experimental uncertainties by the 
ATLAS data 237 . The data are confronted with various pQCD-based model calculations. 


The state-of-the are FONLL pQCD prediction |83|, shown as a solid band in Fig. 33, agrees 


well with the data within experimental and theoretical uncertainties. In order to simplify the 
quantitative comparison of model predictions with the data the lower panels in Fig. [33] show the 
ratio of the data and the other model calculations to the FONLL prediction on a linear scale. 


The NLO calculation POWHEG 238,239 stops on the parton level and, therefore, has to be 


interfaced with parton shower generators to allow a comparison with decay lepton data. Here, 
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Figure 34: p t - and rapidity-differential production cross sections of muons from heavy-flavor 
hadron decays measured with the ALICE experiment in the ranges 2 < p t < 12 GeV/c and 
—4 < y < —2.5 in pp collisions at y/i = 7 TeV (left and middle panel, reprinted from Ref. j22l]), 
and /^-differential cross section at \fs = 2.76 TeV (right panel, reprinted with permission from 
Ref. 1222]; Copyright (2012) by the American Physical Society) in comparison with predictions 
from FONLL pQCD calculations. Error bars (boxes) represent the statistical (systematic) 
uncertainties, where an additional 3.5% absolute normalization uncertainty is not shown. Grey 
bands indicate the FONLL prediction 83 for muons from heavy-flavor hadron decays and 


the long dashed (short dashed) curves represent the FONLL predictions for charm (bottom) 
hadron decays separately. The ratios of the data to the FONLL prediction are shown in the 
lower panels on a linear scale. 


PYTHIA 11201 and HERWIG 12401 have been used for the parton shower simulation. Whereas 
POWHEG+PYTHIA agree reasonably well with the data and with FONLL, this is different 
for POWHEG+HERWIG which predicts a significantly reduced lepton cross section. It was 
shown that this discrepancy can not be accounted for by the different hadron decay models. In 
addition, a LO PYTHIA pQCD calculation 1120 , which includes the parton shower simulation, 
predicts total cross sections approximately a factor two larger than the measurement, but with 
the proper shape of the decay lepton p t distribution. 

It is an interesting question whether the next-to-leading-log (NLL) resummation terms 
included in the FONLL pQCD calculation are actually important at the LHC or whether an 
NLO calculation without this resummation would be sufficient for the description of open heavy- 
flavor data. In order to address this issue, an NLO central value prediction was calculated with 
the FONLL code excluding the NLL resummation |83j. An indication that this approach yields 
too hard lepton p t spectra is visible in the left panel of Fig. |33| already, which shows the electron 
and muon p t spectra up to 26 GeV/c. The right panel of Fig. [33] displays the /^-differential 
production cross section of muons from heavy-flavor hadron decays in the extended transverse 
momentum range up to 100 GeV/c [237|. While the FONLL pQCD prediction and, to some 


extent, also the other pQCD calculations are in good agreement with the measurement, the 
NLO FONLL calculation ignoring the NLL resummation is clearly not consistent with data, 
demonstrating for the first time at a hadron collider the importance of the resummation of high 
order terms in the pQCD expansion. 

At forward rapidity (—4 < y < —2.5) the p t - and rapidity-differential production cross 
sections of muons from heavy-flavor hadron decays have been measured in pp collisions at Qi = 
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7 TeV 2211 and 2.76 TeV |222| with the ALICE experiment as shown in Fig. 34 Corresponding 
FONLL pQCD calculations are in good agreement with the measurement within the systematic 
uncertainties. In both semileptonic decay channels, i.e. electrons and muons, the FONLL pQCD 
calculations suggest that for p t > 5 — 6 GeV/c the dominant contribution is from bottom hadron 
decays while at lower p t charm hadron decays represent the relevant lepton source. 

The correlation of single liigli-/> t . muons with jets was measured with the CMS experiment 

241 . Events with liigli-/> t . muons (+ > 6 GeV/c) have 


in pp collisions at sfs = 7 TeV 
been selected in the range |+j < 2.1. Events containing bottom hadrons were discriminated 
from background via the measurement of the muon transverse momentum with respect to 
the closest jet. The resulting p t - and pseudorapidity-differential cross sections of muons from 
bottom-hadron decays have been confronted with pQCD calculations executed in the MC@NLO 
framework [230] 


While the measured cross section are slightly larger than theoretically pre¬ 


dicted, the data still agree with the model within the combined systematic uncertainties. 


4.2.4 FONLL-driven scaling 

The heavy-flavor production cross sections measured in pp collisions at = 7 TeV can be 
scaled to = 2.76 TeV to provide a reference for the Pb+Pb data at the same energy per 
nucleon-nucleon pair for cases in which a direct measurement of a pp reference at 2.76 TeV is 
not available. This was demonstrated for the ALICE case in Ref. |242|. Since FONLL pQCD 
calculations 831 have been shown to be in reasonable agreement with all open heavy-flavor 


observables measured in pp collisions at \fs = 7 TeV, FONLL was used for the necessary +7 
scaling in this case. The scaling factors for D mesons, electrons, and muons were defined as 
the ratios of the corresponding cross sections from FONLL calculations at 2.76 and 7 TeV. 
The assumption was used that neither the factorization and renormalization scales nor the 
heavy-quark masses in the FONLL calculation vary with \fs. To evaluate the uncertainties 
of the scaling factors, the scales and heavy-quark masses were varied over reasonable ranges 
and the envelopes of the resulting scaling factors were determined. For all open heavy-flavor 
observables the resulting uncertainties of the scaling functions are similar. They range from 
~ 25% at p t = 2 GeV/c to less than 10% for p t > 10 GeV/c. 

As mentioned above, with all LHC experiments pp data have been recorded at = 
2.76 TeV, albeit with limited statistics only. p t -differential muon spectra from heavy-flavor 
hadron decays have been measured with the ALICE experiment 12221 as shown in the right 
panel of Fig. [34} Furthermore, D-meson p tj - differential production cross sections have been 


measured with ALICE at the same energy 1243). Both the muon and the D-meson data are 


in good agreement within uncertainties with the reference obtained via the FONLL-driven \/s 
scaling procedure from the 7 TeV pp measurements to 2.76 TeV. 


4.3 Heavy flavor in Pb+Pb collisions 

4.3.1 Hadronic heavy-flavor hadron decays 

First results on open heavy-flavor production in Pb+Pb collisions at Q+n = 2.76 TeV at the 
LHC have been published recently. One of the major deficiencies of the measurements at RHIC 
has been overcome at the LHC right away. With ALICE, D mesons could be reconstructed even 
in central Pb+Pb collisions 1244] due to the excellent performance of the apparatus. Charged 
particle tracks are reconstructed in the ALICE TPC and ITS with a typical efficiency of 90%. 
Furthermore, the ITS provides excellent vertexing precision, e.g. the track impact parameter 
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resolution in the transverse plane is better than 50 /im for p t > 2 GeV/c. With the TOF system, 
which is instrumental for hadron identification, a time resolution in the order of 100 ps can be 
achieved. The D-meson transverse momentum spectra measured with ALICE in pp and Pb+Pb 
collisions at the LHC contain feed down from bottom-hadron decays. This contribution was 
estimated based on the bottom production cross sections predicted from the FONLL calculation, 
which in pp collisions is in good agreement with the ALICE measurements of J/^ and electrons 

The feed down 


from bottom-hadron decays as demonstrated in Figs. [31| and |32| respectively, 
from bottom hadron decays was subtracted from the inclusive D-meson spectra resulting in 
the so-called prompt D-meson yields. For Pb+Pb collisions, the unknown nuclear modification 
factor of bottom hadrons introduces an additional systematic uncertainty on the prompt D- 
meson yields in Pb+Pb collisions. The conservative assumption was employed that the nuclear 
modification factor of bottom hadrons relative to the prompt D-meson Raa fulfills the condition 
0.3 < Raa/Raa < 3. The resulting nuclear modification factors of prompt D°, D + , and D+ 
mesons in the 20% most central Pb+Pb collisions are shown in the upper left panel of Fig. [35| 
A strong suppression is observed for all species reaching a factor 3-5 for p t > 8 GeV/c. Towards 
more peripheral Pb+Pb collisions nuclear modification factors closer to one are measured for 
prompt D mesons as demonstrated for 40-80% central Pb+Pb collisions in the upper right panel 
of Fig. 35} The centrality dependence of the D-meson Raa is quantified in the middle panel of 
Fig. 35, which shows the nuclear modification factor of D° mesons in the transverse momentum 


range 2 < p t < 5 GeV/c as a function of the average number of participants (N part ) on the left 
and the same dependence for all D-meson species in the range 6 < p t < 12 GeV/c on the right. 
Clearly the D-meson suppression increases towards more central Pb+Pb collisions. 

It is an important issue whether the observed suppression of prompt D mesons at inter¬ 
mediate and high p t might be related to shadowing of the parton distribution functions at 
the LHC. The comparison of an NLO pQCD calculation including the state of the art EPS09 
parametrization of parton shadowing 12451 with the measured D-meson Raa as shown in the 
lower left panel of Fig. [35] provides evidence that this is not the case, i.e. the observed sup¬ 
pression clearly is a final state effect related to the hot and dense medium produced in Pb+Pb 
collisions at the LHC. While parton shadowing seems not to be relevant for the interpretation 
of the prompt D-meson Raa, it is obvious that a highly desirable experimental confirmation of 
this conclusion will require data from proton-nucleus collisions at the LHC. 

One of the most important experimental question left open by heavy-flavor measurements 
at RHIC is related to the expected mass hierarchy of hadron suppression. Radiative in-medium 
energy loss should be larger for light quark flavors compared to charm and, in particular, 
bottom quarks. However, the measured nuclear modification factor of electrons from heavy- 
flavor hadron decays at high p t was comparable in magnitude to the light-flavor dominated 
neutral pion Raa in Au+Au collisions at RHIC, giving no hint of a quark mass hierarchy. The 
more exclusive measurements from Pb+Pb collisions at the LHC might provide new insight into 
this issue. The nuclear modification factor of prompt D mesons measured 1244] as a function 
of p t with ALICE in the 20% most central Pb+Pb collisions at mid-rapidity at the LHC is 
compared with corresponding ALICE measurements of charged hadrons [246 and non-prompt 
J/%> mesons from bottom-hadron decays as measured with the CMS experiment |247| in the 
lower right panel of Fig. [35} While with the current statistical and systematic uncertainties a 
definite conclusion can not be drawn, the ALICE data hint towards a stronger suppression of 
charged hadrons relative to D mesons, as one would expect from a scenario in which induced 
gluon radiation is the dominant energy loss mechanism for partons propagating through a 
dense and color-charged medium. The nuclear modification factor of non-prompt J/% mesons 
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Figure 35: Nuclear modification factor Raa of prompt D°, D + , and D+ mesons measured with 
the ALICE experiment in the 20% most central (upper left) and in 40-80% central (upper right) 
Pb+Pb collisions at ^%nn = 2.76 TeV. Error bars (boxes) show the statistical (systematic) 
uncertainties. The normalization uncertainties are indicated by the full boxes at one. Centrality 
dependence of Raa for D° mesons with 2 < p t < 5 GeV/c (middle left) and all D-meson species 
with 6 < p t < 12 GeV/c (middle right). Average D-meson Raa in the 20% most central Pb+Pb 
collisions in comparison with results from an NLO pQCD calculation including EPS09 nuclear 
shadowing j2451 (lower left). D-meson Raa compared with corresponding values of the nuclear 
modification factors of charged hadrons as measured with ALICE 1246 j and of non-prompt J/+ 
mesons from bottom-hadron decays as measured with the CMS experiment [247] (lower right) 
(reprinted from Ref. 1244] with kind permission from Springer Science and Business Media). 
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Figure 36: Average nuclear modification factor Raa of various D-meson species measured as 
a function of p t with the ALICE experiment at mid-rapidity in the 20% most central Pb+Pb 
collisions at %.snn = 2.76 TeV |244| (left panel) and the corresponding Raa of charged hadrons 
in a slightly wider pseudoradidity interval 1246] (right panel). Error bars (boxes) correspond to 
the statistical (systematic) uncertainties of the measurements. The normalization uncertainties, 
which are almost fully correlated for the two measurements, are indicated by the boxes around 
one. Predictions from various model calculations are compared to the data (see text for further 
details) (reprinted from Ref. 


244 with kind permission from Springer Science and Business 


Media). 


from bottom-hadron decays measured with the CMS experiment, clearly supports this picture 
because it is larger than the Raa measured for charged hadrons. With CMS the Raa of non¬ 
prompt J/xp mesons from bottom-hadron decays was also measured in the 80% most peripheral 
collisions at %.S'nn = 2.76 TeV 247 . Within uncertainties the measured value is consistent 
with the measurement in central Pb+Pb collisions. In addition to the obvious necessity to 
collect statistically more significant data samples it is also mandatory to measure the nuclear 
modification factor of heavy- and light-flavor observables in proton-nucleus collisions at the 
LHC, such that cold nuclear matter effects can be separated from medium modifications due 
to the interaction of heavy quarks with the hot QCD matter produced in nuclear collisions at 
the LHC. 

In the same context, it is important to confront theoretical predictions with the nuclear 
modification factors as measured in central Pb+Pb collisions at the LHC. A comparison of 
various model predictions with the average D meson Raa 244] (left panel) and with the Raa 
of charged hadrons 12461 (right panel) measured with the ALICE experiment in the 20% most 
central Pb+Pb collisions is shown in Fig. 36, In several models the nuclear modification factors 
are calculated both for D mesons and for charged hadrons. Among those, calculations which 
supplement radiative energy loss with the in-medium dissociation of D meson (curves I [I03p48l 
in Fig. [36]) or with collisional energy loss via the WHDG (curves II |249|) or CUJET1.0 (curves 
VII 12501) implementations describe Raa of D mesons and hadrons reasonably well at the 
same time. In the model incorporating in-medium formation and dissociation of D mesons the 
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Figure 37: Elliptic-flow strength v 2 of D° mesons measured with the ALICE experiment 256 
as a function of p t in mid-central Pb+Pb collisions in comparison with the corresponding v 2 of 
unidentified charged particles |65) (reprinted from Ref. [256)). 


medium density is tuned such that it describes the inclusive jet suppression measured at the 
LHC. This is different for the WHDG and CUJET1.0 calculations which extrapolate the medium 
density required to describe the pion suppression observed at RHIC to LHC conditions. This 
extrapolation might contribute to the fact that the nuclear modification factors calculated for 
charged hadrons in these two models are slightly smaller than measured. The average D-meson 
Raa calculated in a completely different model approach based on AdS/CFT drag coefficient 
(curves III [l04j) underestimates the measurement significantly. The predictive power for the 
charged hadron Raa is only limited in this model. Other models, including a Langevin transport 
calculation (curve IV 251,252|), a calculation including collisional and radiative energy loss in 
the medium (curve V |196,253)), a pQCD based partonic transport calculation (curve VI 12541), 


and a calculation employing radiative energy loss only (curve VIII |255]), provide the average 
D-meson Raa only. In general, these model calculations are in reasonable agreement with the 
data with a tendency to underpredict the measured D-meson Raa- 

The measurement of nuclear modification factors of D mesons demonstrates the strong 
interaction of heavy quarks produced in the earliest phase of Pb+Pb collisions at the LHC 
with the hot and dense partonic medium that is formed afterwards. It is an important question 
whether heavy quarks participate in the collective expansion of this medium and whether they 
might thermalize. This question can be addressed via the measurement of the elliptic flow 
strength v 2 of particles carrying heavy quarks. v 2 was measured for neutral D mesons 12561 in 
mid-central Pb+Pb collisions with the ALICE experiment as shown in Fig. [37) A correction 
for feed down from B-meson decays was not applied yet. While statistical and systematic 
uncertainties are still large, an indication for a non-zero v 2 of D mesons is observed. Clearly, 
more statistics is necessary to be able to evaluate the centrality dependence of D-meson flow, 
and it will be important to reduce the systematic uncertainties before any strong conclusion 
can be drawn. 


4.3.2 Semileptonic heavy-flavor hadron decays 

Both electrons and muons from heavy-flavor hadron decays have been measured with the ALICE 
experiment in Pb+Pb collisions at ^/snn = 2.76 TeV. The electron analysis in the p t range up 
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Figure 38: Nuclear modification factor of muons from heavy-flavor hadron decays measured 
as a function of p t with the ALICE experiment in the rapidity range 2.5 < y < 4.0 in the 
10% most central and in 40-80% central Pb+Pb collisions (left panel). Centrality dependence 
of the muon Raa as a function of the average number of participating nucleons in the ranges 
2.5 < y < 4.0 and 4 < p t < 10 GeV/c (right panel) (reprinted with permission from Ref. 1222 ; 
Copyright (2012) by the American Physical Society). 


to 6gevc at mid-rapidity (|y < 0.5) is currently being finalized. Preliminary results 257 


indicate a strong suppression in central collisions consistent with the suppression observed for 
muons from heavy-flavor hadron decays (see below) although the systematic uncertainties are 
still large. Towards more peripheral collisions the high-p t heavy-flavor electron suppression 
decreases, again consistent with the muon case. 

For muons from heavy-flavor hadron decays the nuclear modification factors as measured 


with the ALICE experiment 222 at forward rapidity (2.5 < y < 4.0) in the 10% most central 
and 40-80% central collisions are presented as a function of p t in the left panel of Fig. [38| In 
contrast to the D-meson Raa measurement the muon pp reference could be taken directly from 
the ALICE pp measurement at y/s = 2.76 TeV and not from the FONLL-driven scaling 1242 


of the pp data at 7 TeV. This was possible only because a hardware muon candidate trigger 
allowed a significant muon sample to be recorded in the relatively short pp run period of the 
LHC at 2.76 TeV. In the transverse momentum range 4 < p t < 10 GeV/c a strong suppression 
of muons from heavy-flavor hadron decays is observed in the most central Pb+Pb collisions 
without a pronounced p t dependence. In 40-80% central Pb+Pb collisions a suppression is 
observed as well, although to a much lesser extent, again without any significant p t dependence. 
The pronounced centrality dependence of the muon Raa from heavy-flavor hadron decays is 
quantified in the right panel of Fig. [38] which presents the nuclear modification factor of such 
muons measured with ALICE in the range 6 < p t < 10 GeV/c at forward rapidity as a function 
of the average number of participating nucleons in various Pb+Pb centrality classes. While in 
peripheral Pb+Pb collisions only a small suppression is observed, it increases strongly towards 
central collisions. However, it is important to note that for the decay muons in the p t range 
investigated here, a significant fraction originates from the decay of bottom hadrons. Therefore, 


these data emphasize again (see also the lower right panel of Fig. 35) that not only charm but 
also bottom quarks are coupled strongly to the hot and dense QCD medium produced in high- 
energy heavy-ion collisions. 
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5 Summary and outlook 


Hadrons carrying heavy flavor, i.e. charm or bottom quarks, have been shown to be among 
the most interesting probes in modern particle and nuclear physics. Because of the large 
quark masses heavy-flavor production proceeds mainly through hard parton-parton scattering 
processes in the earliest stage of hadronic collisions. As such the measurement of heavy-flavor 
production in pp collisions provides a crucial testing ground for perturbative QCD calculations. 
State of the art FONLL pQCD calculations are in agreement with all of such measurements 
within statistical and systematic uncertainties at the Fermilab Tevatron, the BNL RHIC, and at 
the CERN LHC. Furthermore, these measurements serve as a baseline for heavy-flavor studies 
in nucleus-nucleus collisions, in which the heavy quarks propagate through and interact with 
the hot and dense medium produced in the nuclear collisions. The investigation of in-medium 
modifications of heavy-flavor observables can shed light on the properties of the QCD medium 
and the nature of parton-medium interactions. This is the case both for heavy-quarkonia, which 
have not been discussed in this review (see Ref. |80| instead) and for which open heavy-flavor 
measurements can provide a natural baseline, and for open heavy-flavor observables, which 
have been in the focus of this review. 

While extensive systematic studies of the production of heavy quarkonia have been con¬ 
ducted in the nucleus-nucleus collision program at the CERN SPS, the production of open 
heavy-flavor hadrons does not play a role in the collision dynamics in this rather low energy 
regime. 

First systematic studies of open heavy-flavor production in nucleus-nucleus collisions have 
been conducted at RHIC, mainly in the semi-electronic decay channel, where it was verified 
that the total heavy-flavor yields scale with the number of binary collisions. Measurements 
of the nuclear modification factor and of the elliptic flow strength of electrons from heavy- 
flavor hadron decays demonstrate the strong coupling of heavy quarks with the hot and dense 
QCD medium produced in Au+Au collisions at RHIC. To further elucidate the properties of 
hot QCD matter it is important to shed light on the relevant interaction mechanism of heavy 
quarks with the medium. It is still not clear which mechanisms are finally responsible for the 
strong suppression observed for high p t electrons in central Au+Au collisions at RHIC. 

At the CERN LHC, the systematic study of open heavy-flavor production in nucleus-nucleus 
collisions is continued in an unprecedented high-energy regime. In addition, for the first time, 
exclusive charm and bottom measurements have become possible in nucleus-nucleus collisions 
thanks to the availability of silicon vertex spectrometers that allow the separation of secondary, 
displaced heavy-flavor hadron decay vertices from the primary collision vertex. The general 
picture established at RHIC was confirmed by the first measurements at the LHC already, 
but major additional progress is expected in the near future through the systematic study of 
the centrality, p t , and rapidity dependence of exclusive charm and bottom observables both 
in Pb+Pb and soon in p+Pb collisions. Table [3] gives an overview over the measured open 
heavy-flavor channels and their kinematic coverage at RHIC and at the LHC until July 2012. 
It should be noted that since then preliminary results have become available, in particular from 
the LHC, which significantly extend the kinematic coverage in most channels. 

The future is bright for open heavy-flavor measurements in high-energy nucleus-nucleus 
collisions. The experimental program at the LHC has just begun. At RHIC, after the very 
successful first decade of operations at the world’s first dedicated heavy-ion collider, the ac¬ 
celerator is upgraded to deliver higher luminosities and concurrently the PHENIX and STAR 
detectors are undergoing major upgrade programs. Key elements of these upgrades are silicon 
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Table 3: Overview over measurements related to open heavy-flavor production conducted at 
RHIC (^/snn < 0.2 TeV) and at the LHC (^/snn < 7 TeV) until July 2012. For each collision 
system, center-of-mass energy ^/snn, and channel investigated the phase space coverage is given 
as well. 


system 

v^nn [TeV] 

channel 

p t coverage [GeV/c 

y coverage ref. 

PP 

0.2 

D° 

0.6 < p t < 2.0 

\y\ < 1 

170 




D* 

2 < p t < 6 

\y\ < 1 

170 




c,b^e ± 

0.3 < p t < 9.0 

\y\ < 0.35 

155 




c,b ^ e ± 

1.2 <p t < 10.0 

\y\ < 1 

174 




c —y e^~ 

2 < p t < 7 

\y\ < 0.35 

163 




b —y 

2 < p t < 7 

\y\ < 0.35 

163 




c —y 

2.5 < p t < 9.5 

12/1 <0.7 

175 




6 —^ e =b 

2.5 < p t < 9.5 

12/1 <0.7 

175 




c, b ^ e + e“ 

Pt > 0.2 

2/ < 0.35 

164 

,165: 



c, b —» 

1 < p t < 7 

1.4 < p < 1.9 

156 


d-H - Au 

0.2 

D° 

0.1 < p t < 3.0 

\y\ < 1 

169 




c,b —>• 

0.85 < p t < 8.5 

2/ < 0.35 

181 




c,b —>■ 

1 < p t < 10 

\y\ < 1 

173 


Cu+Cu 

0.2 

c,b —>• 

1 < Pt < 4 

1.4 < \y\ < 1.9 

156 


Au+Au 

0.13 

c,b ^ e ± 

0.5 < p t < 3.0 

2/ < 0.35 

157 


Au+Au 

0.2 

D° 

Pt < 2.0 

M < i 

171 




c,b e ± 

0.3 < p t < 9.0 

2/ < 0.35 

155 




c,b —>• 

1.2 <p t < 10.0 

M < i 

173 




c,b —>• e + e 

Pt > 0.2 

2/ < 0.35 

165 


PP 

7.0 

D°, D+, D* + 

D + , D* + : 1 < p t < 24 

12/1 <0.5 

2l0 





D°: 1 < p t < 16 





B°, B+, B° 

B+, B°: p t > 5 

B+, B°: |p| < 2.4 

227 

-229 




B°: p t > 5 

B°: \y\ < 2.2 




B+ 

0 < p t < 40 

2.0 <y< 4.5 

231 




b hadron —>■ J /ij> 

Pt > o 

12/1 <0.9 

211 




b hadron —> 3/xp 

7 < pt < 70 

2/ < 0.75 

234 




b hadron —> 3/xp 

1 < pt < 30 

1.5 < |y| < 2.0 

234 




b hadron —> 3/ip 

Pt > 6.5 

\y\ < i-2 

235 




b hadron —> 3/ip 

Pt > o 

1.6 < \y\ < 2.4 

235 




c,b —>• 

0.5 < pt < 8.0 

12/1 <0.5 

212 




c,b ^ e ± 

7 < pt < 26 

\y\ < 2.0 

237 




c —» e ± 

1 < pt < 8 

12/1 <0.8 

220 




b —y 

1 < p t < 8 

12/1 <0.8 

220 




c,b —)■ 

2 < p t < 12 

2.5 < y < 4.0 

221 




c,b —>• 

7 < p t < 24 

2/ < 2.0 

237 


PP 

2.76 

D°, D+, D* + 

1 < p t < 12 

12/1 <0.5 

243 




c, b —^ 

2 < p t < 10 

2.5 < y < 4.0 

222 


Pb+Pb 

2.76 

D°, D+, D* + 

2 < p t < 16 

12/1 <0.5 

243 




b hadron —> J/-0 

6.5 < pt < 30 

2/ < 2.4 

247 




c, 6 —» e 1 * 1 

3.5 < pt < 6.0 

12/1 <0.5 

257 




c, b —^ 

4 < p t < 10 

2.5 < |y| < 4.0 

222 



Ttr 



























































































vertex detectors aimed towards the separation of displaced heavy-flavor hadron decay vertices 
from the primary collision vertex. With these upgrades in place it will be possible both at 
RHIC [152 as well as at the LHC to measure the yields and phase space distributions of charm 
and bottom hadrons in pp, p/d-nucleus, and nucleus-nucleus collisions via exclusive channels, 
providing unique hard probes for hot QCD matter even at low p t . In particular the measure¬ 
ment of heavy-flavor baryons, notably the A c , would be important even though this is very 
challenging, markedly in nucleus-nucleus collisions. 

Precision measurements of the nuclear modification factors of heavy-flavor hadrons will 
address the issue whether the expected mass ordering of quark energy loss in the hot and 
dense medium is realized in nature or not. This is one of the crucial questions related to the 
in-medium parton dynamics and, in particular, to the mechanism of parton energy loss in the 
medium. The latter is up to now not understood on a fundamental level. Tightly connected to 
the energy loss of a hard probe is the reaction of the medium to the resulting energy deposit. 
Regarding these issues, heavy-flavor measurements can contribute a unique piece of information. 
With the high resolution vertex detectors at RHIC and at the LHC it should be possible to 
“tag” charm and, in particular, bottom jets on an event-by-event basis. In contrast to high-p t 
particle production, which at RHIC and at the LHC is dominated by gluon jets, such heavy- 
flavor tags give clean access to quark jets. Furthermore, momentum conservation in the hard 
scattering process gives rise to a back-to-back correlation in azimuth of the quark and antiquark 
jets produced associately. Therefore, tagged heavy-quark jets constitute a new quality of hard 
probes as they give much better access to the kinematics of the initially scattered partons and 
allow the investigation of the mass dependence of parton energy loss and the medium response. 
Also, the precision measurement of heavy-flavor hadron flow will be possible at RHIC and at the 
LHC, such that the question can be addressed to what extent heavy quarks reach equilibrum 
with the hot and dense medium. This is a crucial issue as it is directly related to the transport 
properties of the medium. 

With precision vertexing it should furthermore be possible to measure at RHIC and at the 
LHC the contribution of correlated semileptonic heavy-flavor hadron decays to the dilepton 
continuum in the intermediate mass region, which at RHIC is dominated by charm hadron 
decays while at the LHC bottom hadron decays contribute significantly. In nucleus-nucleus 
collisions, these dileptons carry information about the interaction of heavy quarks with the 
medium in a twofold way. First, the quarks lose energy in the medium and second the opening 
angle between the quark and the antiquark can be modified. In the extreme case, the heavy 
quarks thermalize with the medium such that the angular correlation between the quark and 
antiquark is completely lost. Again, a comparison of charm and bottom measurements would 
be useful regarding the issue of the parton energy loss mechanism and it could help to better 
understand the transport properties of the medium. Furthermore, it should be noted that 
only two sources are expected to contribute significantly to the intermediate mass dilepton 
continuum, i.e. correlated heavy-flavor hadron decays and thermal radiation from the hot 
medium. If these two contributions can be disentangled with the help of precision vertexing 
the temperature of the emitting source and its evolution with time can be investigated in a 
unique way. 


In the longer term future, open charm production will be addressed at the FAIR facility 
currently under construction at GSI. At this fixed-target machine charm can be investigated 
close to the production threshold in proton-nucleus and nucleus-nucleus collisions, where models 
predict in-medium effects which in a unique way could shed light on the properties of the 
baryon-rich form of dense QCD matter produced at FAIR 1129 . 
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